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Abstract 

Wc reexamine and further develop different gravito-electromagnetic (GEM) analogies found 
in the literature, and clarify the connection between them. Special emphasis is placed in two 
exact physical analogies: the analogy based on inertial fields from the so-called "1-1-3 formalism", 
and the analogy based on tidal tensors. Both are reformulated, extended and generalized. We 
write in both formalisms the Maxwell and the full exact Einstein field equations with sources, 
plus the algebraic Bianchi identities, which are cast as the source-free equations for the gravita- 
tional field. New results within each approach are unveiled. The well known analogy between 
linearized gravity and clcctromagnetism in Lorentz frames is obtained as a limiting case of the 
exact ones. The formal analogies between the Maxwell and Weyl tensors, and the related issue 
of super-energy, are also discussed, and the physical insight from the tidal tensor formalism 
is seen to yield a suggestive interpretation of the phenomenon of gravitational radiation. The 
precise conditions under which a similarity between gravity and electromagnetism occurs are 
discussed, and wc conclude by summarizing the main outcome of each approach. 
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1 Introduction 

This work has two main goals: one is to establish the connection between the several gravito- 
electromagnetic analogies existing in the literature, summarizing the main results and insights 
brought by each of them; the other is to further develop and extend some of these analogies. 

In an earlier work by one of the authors [H [2] , a gravito-electromagnetic analogy based on tidal 
tensors was presented, and its relationship with 1) the well known analogy between linearized gravity 
and electromagnetism, 2) the mapping, via the Klein-Gordon equation, between ultrastationary 
spacetimes and magnetic fields in curved spacetimes, and 3) the formal analogies between the 
Weyl and Maxwell tensors (their decomposition in electric and magnetic parts, the quadratic scalar 
invariants they form, and the field equations they obey) was discussed. 

Building up on the work in [2], another approach is herein added to the discussion: the exact 
analogy based on the fields of inertial forces, arising in the context of the 1+3 splitting of spacetime. 
This approach, which is herein reformulated and suitably generalized, is still not very well known. 
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but very far reaching. It is therefore important to understand how it relates with the other known 
analogies, and in particular with the (also exact) approach based on tidal tensors. 

Each of the analogies discussed here are also further developed, and some new results within 
each of them are presented. We start in Sec. [2] by revisiting the analogy based on tidal tensors 
introduced in [1] (and partly reviewed in [3]). In [1] it was shown that there is an exact, physical 
gravito-electromagnetic analogy relating the so-called electric (E^^) and magnetic (Hq,^) parts of 
the Riemann tensor with the electromagnetic tidal tensors {Eai3, Bap) defined therein: the "relative 
acceleration" between two nearby test particles, with the same ^-velocity in a gravitational field, 
is given by a contraction of E^^^ with the separation vector 5x^; just like the "relative acceleration" 
between two nearby charged particles (with the same U'^ and ratio q/m) is given by a contraction 
of the electric tidal tensor E^p with fe". Moreover, the gravitational force exerted on a spinning 
particle is exactly given by a contraction of Hq,^ with the spin vector 5" (the "gravitomagnetic 
dipole moment"), just like its electromagnetic counterpart is exactly given by a contraction of the 
magnetic tidal tensor Bap with the magnetic dipole moment Here Eq,^, EIq,/3, i?Q,^ are the 
tidal tensors as measured by the test particles. In this work we unveil another exact physical analogy 
involving Ba/s and Ma/s, manifest in the equations for the relative precession of spinning particles: 
both in electromagnetism and gravity, the relative precession of a test gyroscope/magnetic dipole 
relative to a system of axes attached to guiding gyroscopes/magnetic dipoles at a neighboring point, 
moving with the same 4-velocity U", is given by a contraction of Ba/s /^ai3 with the separation vector 
6x^ . This gives a physical interpretation of the tensors Ba/^ and complementary to the one 
described above. The expression for the differential precession of gyroscopes in terms of M^js was 
originally found in a recent work |102j : herein (Sec. 12.3]) we re-derive this result through a different 
procedure (which, we believe, is more clear), and we obtain its electromagnetic analogue. 

The analogy based on tidal tensors extends to the field equations of both theories; in [1] it was 
shown that by taking the traces and antisymmetric parts of the electromagnetic tidal tensors one 
obtains Maxwell equations, and performing the same operations in the gravitational tidal tensors 
leads to a strikingly similar set of equations, which turn out to be the temporal part of Einstein's 
equations. Herein, building on that work, we extend this formalism to the full gravitational field 
equations (Einstein equations with sources plus the algebraic Bianchi identities); we show that 
they can be decomposed in a set of equations involving only the sources, the gravitoelectric E^^ 
and gravitomagnetic Hq,^ tidal tensors, plus a third spatial tensor Fq^, introduced by Bel |48| . 
which has no electromagnetic analogue. More precisely, making a full 1-1-3 covariant splitting of 
the gravitational field equations, one obtains a subset of four equations (involving the temporal 
part of the curvature), which are formally similar to Maxwell equations written in this formalism 
(as found in [1]), both being algebraic equations involving only tidal tensors and sources, plus an 
additional pair of equations involving the purely spatial curvature (encoded in Eq^), which have no 
electromagnetic counterpart. In the case of vacuum, this latter part identically vanishes, and not 
only the full Einstein equations become algebraic equations for tidal tensors, just like the Maxwell's, 
but also one can actually exactly obtain the former by simply replacing {Eajj, Bap} — )• {Eo,^,EIq,^} 
in the latter. It is discussed how this approach, by bringing together the two theories in a single 
formalism, is especially suited for a transparent comparison between the two interactions, which 
reduces to comparing the tidal tensors of both sides (which is straightforward in this framework). 
Fundamental differences between the two interactions are encoded in the symmetries and time 
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projections of the tidal tensors (as already pointed out in [T]); herein we explore the consequences 
in terms of the worldline deviation of (monopole) test particles; and in the companion paper [3], in 
the dynamics of spinning multipole test particles. 

In Sec. [3l we discuss another exact gravito-electromagnetic analogy, the one drawing a paral- 
lelism between spatial inertial forces (described by the GEM fields) and the electromagnetic fields. 
GEM fields are best known from linearized theory, e.g. |23l \98\ [971 124] : less well known is the 
exact "Quasi-Maxwell" analogy, based on the inertial fields that arise in the 1+3 splitting of sta- 
tionary spacetimes introduced by Landau-Lifshitz [77j . and further worked out by many authors 
[TH| [HI [ini [Ml [Ml [M] • Even less known is the existence of an exact formulation applying to arbi- 
trary spacetimes [191 120j. which considers an arbitrary timelike congruence and does a general 1+3 
splitting in terms of a time direction parallel to the corresponding observers' worldlines, and the 
local rest spaces orthogonal to them. In this latter framework, the gravito-electromagnetic analogy 
(in the sense of a one to one correspondence) is, for most effects, only recovered for stationary 
spacetimes (exceptions are the case of the spin precession, and the hidden momentum of a spinning 
particle discussed in [3]). Herein we reformulate this approach in a slightly more general form, in 
the sense that we use an arbitrary (not the congruence adapted) frame, which splits the gravito- 
magnetic field H in its two constituent parts, of different mathematical origin: the vorticity w of 
the observer congruence and the rotation 0, of the local tetrads (associated to each local observer) 
relative to local Fermi- Walker transport. This degree of generality is of use both in the compan- 
ion paper [5], and in [M]- It also allows, for instance, one to realize that there is a distinction 
between the gravitomagnetic effects detected in the analysis of the LAGEOS Satellites data [41j 
(and presently under scrutiny in the ongoing LARES mission |43] ) . and the one measured by the 
Gravity Probe B experiment. An important result of this approach is an exact expression for the 
geodesic equation, written in the language of GEM fields, that is fully general (valid for arbitrary 
fields, and formulated in terms of an arbitrary frame); the approximate GEM expressions from 
the Post-Newtonian and linearized theory in the literature, as well as the exact expressions of the 
Quasi-Maxwell formalism for stationary spacetimes, are just special cases of this general equation. 
The exact, mathematically rigorous formulation yields a powerful formalism with a very broad spec- 
trum of applications; and an adequate account of some subtleties involved which are overlooked 
in the more common linear approaches. Indeed, the inertial GEM fields in this formulation — 
the gravitoelectric field G, which is but minus the acceleration of the congruence's observers, and 
the gravitomagnetic field H = uj + Q — can be regarded as the general, exact form of the GEM 
3-vectors fields of the usual linearized theory, which (in the way they a are usually presented) are 
somewhat naively defined from the temporal components of the metric tensor (drawing a paral- 
lelism with the electromagnetic potentials), without making transparent their status as artifacts of 
the reference frame, and in particular their relation with the kinematical quantities associated to 
the observer's congruence. Problems which in the approximate descriptions usually end up being 
treated superficially, overlooking the complicated underlying problem, are the effects concerning 
gyroscope "precession" relative to the "distant stars" (such as the Lense-Thirring and the geodetic 
precessions) — the question arising: how can one talk about the "precession" of a local gyroscope 
relative to the distant stars, as it amounts to comparing systems of vectors at different points 
in a curved spacetime? Such comparison can be done in a certain class of spacetimes, and the 
mathematical basis for it is discussed in Sees. [STTI and [331 
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We clarify the relationship between these inertial fields and the tidal tensors of Sec. [2] (a partic- 
ularly important result in the context of this work). For stationary spacetimes, and frames adapted 
to stationary observers (Quasi-Maxwell regime), we express, as in Sec. [21 the Maxwell equations, 
and the full Einstein equations and algebraic Bianchi identities, in this formalism. Again, a set of 
four equations — the time-time and time-space components of the gravitational equations - are 
produced which exhibit a striking formal similarity with their electromagnetic counterparts; and 
there is a fifth equation, the space-space component of Einstein's equations, which has no grav- 
itational analogue. We show also that an exact GEM analogy relating the force on a gyroscope 
with the force on a magnetic dipole exists also in this formalism. In an earlier work by one of the 
authors [H], the force exerted on a gyroscope, at rest with respect to a rigid frame in a stationary 
spacetime was written in terms of the GEM fields associated to that frame; herein we obtain the 
electromagnetic counterpart: the force exerted on a magnetic dipole at rest in a rigid, arbitrarily 
accelerated frame, in terms of the electromagnetic fields associated to that frame, and show there 
is a one to one correspondence with the gravitational analogue. 

In Sec. m we discuss a special class of spacetimes admitting global rigid geodesic congruences, 
the "ultra-stationary" spacetimes. They have interesting properties in the context of GEM, which 
were discussed in [T]: they are exactly mapped, via the Klein-Gordon equation, to magnetic fields 
in curved spacetimes, and have a linear gravitomagnetic tidal tensor, matching the electromag- 
netic analogue. Herein we revisit those spacetimes in the framework of the GEM inertial fields 
of Sec. [31 which sheds light on their properties: these metrics are characterized by a vanishing 
gravitoelectric field and a non- vanishing gravitomagnetic field which is linear in the metric tensor; 
these two properties together explain the above mentioned exact mapping, and the linearity of the 
gravitomagnetic tidal tensor. We also interpret the non-vanishing of gravitoelectric tidal tensor, a 
question left open in [1]. 

In Sec. [5] we explain the relation between the exact approach based in the inertial GEM fields 
of Sec. [31 and the popular gravito-electromagnetic analogy based on linearized theory, e.g. [231 [Ml 
[64l I97j . The latter is obtained as a special limit of the exact equations of the former. Taking 
this route allows, as explained above, to have a clearer account of the physical meaning of the 
GEM fields and other kinematical quantities involved. It is also a procedure for obtaining the 
field equations in terms of (physically meaningful) GEM fields that does not rely on choosing the 
harmonic gauge condition and its inherent subtleties (that have been posing some difficulties in 
the literature, see e.g. [971 11191 [Ml [2])- The usual expression for the force on a gyroscope in the 
literature (e.g. [3 [211 [23]) is also seen to be a (very) limiting case of the exact equation equation 
given in the tidal tensor formalism of Sec. [21 

In Sec.[6]we discuss the analogy between the electric (Safs) and magnetic (Tiafs) parts of the Weyl 
tensor, and the electromagnetic fields E"' and B°'. As already discussed in [2], this is a purely formal 
analogy, as it deals with tidal tensors in the gravitational side, and vector fields (not tidal tensors, 
which are one order higher in differentiation) of electromagnetism. One of the major results of this 
approach is that the differential Bianchi identities (the higher order field equations), when expressed 
in terms of £ai3 and T-Lai3, exhibit some formal similarities with the Maxwell equations written 
in terms of E°', B"', measured with respect to an arbitrarily shearing, rotating and expanding 
congruence. In the case of vacuum, in the linear regime, they become Matte's equations — a set of 
equations formally similar to Maxwell's equations in a Lorentz frame, only with the gravitational 
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tidal tensors {£ij = E^-, Tiij = Mij) in place of the electromagnetic fields. We arrive at a suggestive 
interpretation of gravitational radiation as a pair of traveling orthogonal tidal tensors, propagating 
by mutually inducing each other (just like E and B in the electromagnetic waves), and whose 
effects on test particles are readily established via the tidal tensor approach of Sec. [2l Such "tidal 
tensor interpretation" of the gravitational waves has not been put forth in the literature, possibly 
for two reasons: one is that in the literature dealing with the higher order field equations the 
physical interpretation of Tia/s appears as an unanswered question (HU [87l [56l [82l [83l [32] ; the other 
is that in the literature dealing with the interaction of gravitational waves with spinning particles 
[72\ \7T\ [75t , neither the waves are written in tidal tensor form (but instead in the more usual 
equations for the propagation of the metric perturbations), nor the force on a spinning particle is 
cast in terms of 7iai3, since the Mathisson-Pirani condition is not employed (the latter is necessary 
in order to make TiajS explicitly appear in the equations, see [1]). This analogy has been used 
to address the problem of the definition of local quantities for the gravitational field analogous to 
the electromagnetic field energy and Poynting vector, which has motivated the definition of the 
Bel- Robinson tensor (the so-called "superenergy" tensor). The viewpoint that gravitational waves 
should be thought of as carrying superenergy, not energy, fits well in the interpretation above, and 
with what it tells us about the interaction of waves with multipole test particles. 

1.1 Notation and conventions 

1. Signature and signs. We use the signature — h ++; ea/Sa--/ = ^—g[a/3j5] denotes the 
Levi-Civita tensor, and we follow the orientation [1230] = 1 (i.e., in flat spacetime 
€1230 = !)• ^ijk = ^ijko is the 3-D alternating tensor. ★ = Hodge dual. 

2. Time and space projectors. (T")"^ = —u°'up, (/i")"^ = u"^up + g'^^ are, respectively, the 
projectors parallel and orthogonal to a unit time- like vector u"; may be interpreted as the 
time and space projectors in the local rest frame of an observer of 4- velocity u". (a) denotes 
the index of a spatially projected tensor: A^"^'^"' = {h^)°^pA^^'" . 

3. pc = —j^Ua and j" are, respectively, the charge density and current 4-vector; p = Ta/^lf^U^ 
and J° = —T'^U^ are the mass/energy density and current (quantities measured by the 
observer of 4- velocity U"); T^p = energy-momentum tensor. 

4. = spin 4-vector; p" = magnetic dipole moment, defined as vectors whose components in 
the particle's proper frame are 5" = {0, S), p = (0, p). 

5. Tensors resulting from a measurement process. (|^")"i--"n denotes the tensor A as measured 
by an observer 0{u) of 4-velocity n". For example, (£^")" = F°'^u^ , {E^)a/3 = Fa-^-pu^ 
and (E")c/3 = Rau/Suu'^u^ denote, respectively, the electric field, electric tidal tensor, and 
gravito-electric tidal tensor as measured by 0{u). Analogous forms apply to their mag- 
netic/gravitomagnetic counterparts. 

For 3- vectors we use notation A{u); for example, E{u) denotes the electric 3- vector field as 
measured by 0{u) (i.e., the space part of (£^")", written in a frame where = 0). Often we 
drop the superscript (e.g. (S^)" = -E"), or the argument of the 3- vector: E{U) = E, when 
the meaning is clear. 
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6. Electromagnetic field. The Maxwell tensor F"^^ and its Hodge dual •kF°'^ = e^y F^^ /2 
decompose in terms of the electric (E'")" = F°jjU^ and magnetic (S")" = -kF'^^u^ fields 
measured by an observer of 4-velocity u'^ as 

F^p = 2u^^{E^)p^+e^p^su'{B^V (a) ^ F^^ = 2u^^{B^) - e^p^^u^E'-)'^ (b) (1) 

2 The gravito-electromagnetic analogy based on tidal tensors 

The rationale behind the tidal tensor gravito-electromagnetic analogy is to make a comparison be- 
tween the two interactions based on physical forces present in both theories. The electromagnetic 
Lorentz force has no physical counterpart in gravity, as monopole point test particles in a grav- 
itational field move along geodesies, without any force being exerted on them. In this sense the 
analogy drawn in Sec. 13. 2^ between Eqs. (|46p - (j45p . is a comparison of a physical electromagnetic 
force to an artifact of the reference frame. Tidal forces, by their turn, are covariantly present in 
both theories, and their mathematical description in terms of objects called "tidal tensors" is the 
basis of this approach. Tidal forces manifest themselves in essentially two basic effects: the rel- 
ative acceleration of two nearby monopole test particles, and in the net force exerted on dipoles. 
These notions of multipole moments arise from a description of the test bodies in terms of the 
fields they would produce. In electromagnetism is the multipole expansion of the current density 
vector j° = (pcj), rigorously established in [9], and well known in textbooks as the moments of 
the charge and current densities. In gravity is are the moments of the energy momentum ten- 
sor Tq^, the the so called [8] "gravitational skeleton", of which only the moments of the current 
density J° = —T^^Ufs have an electromagnetic counterpart. Monopole particles in the context of 
electromagnetism are those whose only non-vanishing moment is the total charge; dipole particles 
are particles with nonvanishing electric and magnetic dipole moments (i.e., respectively the dipole 
moments of pc and j); see [9] and companion paper [4] for precise definitions of these moments. 
Monopole particles in gravity are particles whose only non- vanishing moment of T°'^ is the mass, 
and correspond to the usual notion of point test particle, that moves along geodesies. There is no 
gravitational analogue to the intrinsic electric dipole, as there are no negative masses; but in the 
multipole scheme, there is an analogue to the magnetic dipole moment, which is the "intrinsic" an- 
gular momentum (i.e. the angular momentum about the particle's center of mass), usually dubbed 
spin vector/tensor. Sometimes it will also be dubbed herein, for obvious reasons, "gravitomagnetic 
dipole moment". A particle possessing only pole-dipole gravitational moments corresponds to the 
notion of an ideal gyroscope. We have thus two physically analogous effects suited to compare grav- 
itational and electromagnetic tidal forces: worldline deviation of nearby monopole test particles, 
and the force exerted on magnetic dipoles/gyroscopes. An exact gravito-electromagnetic analogy, 
summarized in Table [U emerges from that comparison. 

Eqs. ([Hi) are the worldline deviations for nearby test particles with the sam^ tangent vector 
(and the same ratio charge/mass in the electromagnetic case), separated by the infinitesimal vector 

^We would like to emphasize this point which, even today, is not clear in the literature. Eqs. ([l]l) apply only 
to the instant where the two particles have the same (or infinitesimally close, in the gravitational case) tangent 
vector. This seems to be well understood in the case of the electromagnetic deviation ([T]la), but is not so with the 
geodesic deviation equation ([T]la); which is presented in some literature as applying to the case where the deviation 
velocity is not zero [281 11201 1121j . there being claims that it is an arbitrary initial condition [28]. That is not the 
case: for both electromagnetism and gravity, in the more general case that the velocity of the two particles is not 
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Table 1: The gravito-electromagnetic analogy based on tidal tensors. 



Elcctromagnetism 




Gravity 




Worldline deviation: 


(Pla) 


Geodesic deviation: 


(Elb) 


Force on magnetic dipole: 




Force on gyroscope: 






(1112a) 




(1112b) 


Differential precession of magnetic dipoles: 




Differential precession of gyroscopes: 






tyjod ) 


AO' — n-R* At^ 




Maxwell Source Equations 




Einstein Equations 




• Time Projection: 

= 4^pe 


(lll3a) 


• Time-Time Projection: 


(lll3b) 


• Space Projection: 


(Pea) 


• Time-Space Projection: 


(1116b) 



• Space-Space Projection: 

No electromagnetic analogue F"j + E"_j — F'^^/i"^ = Stt 



Bianchi Identity 




Algebraic Bianchi Identity 




★F"^^^ = ^^[,^.^] = ) 




=0 R[al3-f]S = 0) 




• Time Projection: 




• Time-Time (or Space-Space) Proj: 






(lll5a) 




dllSb) 


• Space Projection: 




• Space-Time Projection: 






(1114a) 


E[a/3] = 


(1114b) 


No electromagnetic analogue 




• Time-Space Projection: 





Jrr". They tell us that the so-called electric part of the Riemann tensor E"^ = R'^^pj^U^U'^ 
plays is geodesic deviation equation dUlb) the same physical role as the tensor Eap = Fa^-pV^ 
in the electromagnetic worldline deviation dUla). E^p measures the tidal effects produced by the 
electric field = FZ^U^ as measured by the test particle of 4- velocity U". We can define it as a 
covariant derivative of the electric field as measured in the inertial frame momentarily comoving 
with the particle: Ea/s = -Eo;/3|c/=const.- Hence we dub it "electric tidal tensor", and its gravitational 

infinitesimally close, the deviation equations include more terms (which depend on both particle's 4-velocity, thus in 
their relative velocity); the relative acceleration is not, in either case, given by a simple contraction of a tidal tensor 
with a separation vector, see 1T]I2]I2H1. A more detailed discussion of this important issue will be presented elsewhere. 
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counterpart "gravitoelectric tidal tensor". Eqs ([Il2) are respectively, the electromagnetic force 
on a magnetic dipole [1], and the Mathisson-Papapetrou-Pirani [HI [lOl [TTl H] equation for the 
gravitational force exerted on a gyroscope. They tell us that the magnetic part of the Riemann 
tensor H"^ = -kR°'^pJJ^U'^ plays in the gravitational force ([U2b) the same physical role as the tensor 
Baj3 = •kFa-y-jsU'^ in the electromagnetic force ([Tj2a). Ba/B measures the tidal effects produced by 
the magnetic field B"' = •kF^U'^ as measured by the particle of 4- velocity for this reason we 
dub it "magnetic tidal tensor", and its gravitational analogue "gravitomagnetic tidal tensor". 

2.1 Tidal tensor formulation of Mcixwell and Einstein equations 

Taking time and space projections Maxwell's and Einstein's can be expressed in tidal tensor formal- 
ism; that makes explicit a striking aspect of the analogy: Maxwell's equations (the source equations 
plus Bianchi identity) may be cast as a set of algebraic equations involving only tidal tensors and 
source terms (the charge current 4- vector); and Einstein's equations (the source equations plus the 
algebraic Bianchi identity) consist of two parts: i) a temporal part that is formally very similar 
to Maxwell's equations and are likewise algebraic equations involving only the gravitoelectric and 
gravitomagnetic tidal tensors and sources (the mass-energy current vector), and ii) a space-space 
part which involves a third spatial tensor and has no electromagnetic analogue. That is what we 
are going to show next. For that we first introduce the time and space projectors with respect to 
a unit time-like vector [/" (i.e., the projectors parallel and orthogonal to 17°'): 

T"^ = {T^r^ = -U'^Up; = {h^r^ = U-Up + 5%. (2) 

A vector A" can be split in its time and space projections with respect to [/"; and arbitrary rank 
n tensor can be completely decomposed taking time and space projections in each of its indices 

(e.g. m)-- 

= T'^pAP + h'^^A^; = (j''^^ + ... (^T"^"^ + (3) 

Instead of using h^^, one can also, if convenient, spatially project an index of a tensor A'^'" con- 
tracting it with the spatial 3- form CaiSa'yU"'', for instance, for the case of vector A'^ , one obtains 
the spatial 2-form Ea^a-yU^ A'^ = •kAajj^U"' , which contains precisely the same information as the 
spatial vector A^h"^^ = A^'^'^ (the former is the spatial dual of the latter). New contraction with 
f^^fivU^ yields A^"^ again. Indeed we may write 

Another very useful relation is the following. The space projection h^ah'^j^F^y = of a 

2-form Fap = is equivalent to the tensor e^yapF'^^U^ = -kF^yU^ (i.e., spatially projecting F^p 
is equivalent to time-projecting its Hodge dual). We have: 

i^H(/3> = e^apxUh^.^sU^'F'^' = e^^p^U^ * F^^.U^ (4) 

Now let -F'yi...'y„o/35i...5^ = -^7i...7„[q/3]5i...<5„; be some tensor antisymmetric in the pair Q,/3; an 
equality similar to the first one above applies: 

-^7i...7„(a>(/3>5i...5^ = ^tio.pxU^e^^^^V F^^^ (5) 

9 



2.1.1 Maxwell's equations 

Maxwell equations are cast in tensor form by the pair of equations: 

^"'';/3 = 4^r (a); *i^"'';;3 = (b). (6) 

being ([H^) the Maxwell source equations and the source free Eq. dSJa), equivalent to F^afS-.^y] = 0) is 
commonly called the Bianchi identity for These equations can be expressed in terms of tidal 
tensors using the decompositions 

Fa/3;'y = 2[/[^£;^]^ + f^aP^iaU" B^^ (7) 
-^Fap-^ = 2U[aBp^^ - eap^aU'^ E^-f (8) 

These expressions are obtained decomposing the tensors Fap-^-y and -kFap-^^ in their time and space 
projections in the first two indices, using Eq. ^ to project spatially. Taking the time projection of 
([6^), we obtain Eq. ([Tj3a) of Table [D taking the space projection by contracting with the spatial 
3-form e^uaaU^ , yields Eq. ([Tj6a). The same procedure applied to Eq. ([6)3), yields Eqs. ([Il4a) and 
([Tj5a) as time and space projections, respectively. 

Hence, in this formalism. Maxwell's equations are cast as the Eqs. of the traces and antisym- 
metric parts of the electromagnetic tidal tensors; and they involve only tidal tensors and sources, 
which is easily seen substituting the decompositions ([7])- ([8]) in Eqs ([Il4a) and ([Il6a), leading to the 
equivalent set: 

E'^a = ^-^Pc (9) 
E[ap] = U^^Ep^^W + ]^e^p^„U''B^^U^ (10) 
B^^ = (11) 
= Uy^Bp^^W -]^e^p^„U''E^'''U.,-2^ie^p,^fW (12) 

The pair of Eqs. ()10p and ()12p can be condensed in the equivalent pair 

^^\sU^E^^p^ = -B^pU^- (a) e^\^U^B^^p^=Eo^pU^ + AT,jo^ (b) (13) 

In a Lorentz frame, since C/"^ = f7°^ = 0, we have E.^^ = E^-js, B^js = B^.p] and (using [7° = 6q) 
Eqs. (|13|) can be written in the familiar vector forms V x E = —dB/dt and V x = dE jdt + 47rj, 
respectively. Likewise, Eqs. ([9]) and pip reduce in this frame to the familiar forms V • -E = 47rpc 
and V • = 0, respectively. 

2.1.2 Einstein's equations 

Equations ()14h ,) below are the Einstein source equations for the gravitational field; Eqs (|14b) are 
the algebraic Bianchi identity, equivalent to R\a^^\& = 0: 

R\-y^^Rai3=8n(^T^l,-^g^pT^^^ (a); = (b). (14) 
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In order to express these equations in the tidal tensor formahsm we will decompose the Riemann 
tensor in its time and space projections (in each of is indices) with respect to a unit time-like vector 
Ija^ j^afS-yS = (ja^ + /^o^) , _ + ^p..a^ gq. ([3]); we obtain 

+6°^*'^t/^6^;,[/,F<^^. (15) 

where we made use of the identity (jl]) to project spatially an antisymmetric pair of indices, noting 
that Rap-yS can be regarded as a double 2-form. This equation tells us that the Riemann tensor 
decomposes, with respect to ?7°, in three spatial tensors: the gravitoelectric tidal tensor Eq,/?, the 
gravitomagnetic tidal tensor M.^^, plus a third tensor 

introduced by Bej§ |48j . which encodes the purely spatial curvature with respect to 11°^, and has 
no electromagnetic analogue. In order to obtain Eq. ()15p . we made use of the symmetries R^l^'^^ = 
RW]b&]^ and in the case of the terms involving Hq/j (and only for these terms) we also assumed 
the pair exchange symmetry R^P'i^ = R'^^^P . 'K^p and Fq,^ are symmetric and therefore have 6 
independent components each; EIc,/3 has 8, together encoding the 20 independent components of the 
Riemann tensor. 

In what follows we will need also the Hodge dual, in the first two indices, of the decomposition 

m- 

^R-^^, = 26"^,^E^[^C/,][/- + 4C/["E[^|,C/,]+e"'\^6^;,[/.]HI^-^^ 

-2U^^¥P\e'"'^^U, . (16) 

The Ricci tensor R^^ = R'^^^g and the tensor ^R^^^^p follow as: 

= -e"^^^M^,[/5C/, - to^s^^WU^U'' - - E''^ + Wf^U^Us + F'^^/i^^, (17) 

* R'^L = e'-^xXaUsU^ - + [/'^F"^6^^„,C/.. (18) 

Substituting (|17p in (|14fe i). and (jlSp in (jl4b ) we obtain Einstein's equations and the algebraic 
Bianchi identities in terms of the tensors Eq,^, Hq^, Fq^. Now let us make the time-space splitting 
of these equations. Eq. (|14b ) is symmetric, hence it only has 3 non-trivial projections: time-time, 
time-space, and space-space. The time-time projection yields 

E-„ = 47r(2/,^ + n) (19) 

where pm = T'^^UpUa denotes the mass-energy density as measured by an observer of 4- velocity 
U^. Contraction of (fT7|) and (fHk ) with the time-space projector 1'^ ^e^ ^^^V^ yields: 

M[,,] = -47reA..^ J^C/^ . (20) 

^The characterization of the Riemann tensor by these three spatial Rank 2 tensors is known as the "Bel decom- 
position"; even though, for some reason, the explicit decomposition (|15l) has not, to the authors' knowledge, been 
written in the literature. 
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The space-space projection yields: 



¥^g + E^g - W^^h^g = 8tt h 




(21) 



where T, ^ ^ h\h%T,'. 



Eq. ()14b ). not exhibiting the tensor -kB?"^^^ a particular symmetry, seemingly splits in four parts: 
a time-time, time-space, space-time, and space-space projections. However only three of them are 
distinct. Substituting decomposition (jlSp in Eq. (jl4b ). and taking the time-time, time-space, and 
space-time projections, yields, respectively: 



The space-space projection yields again H"^, i.e., the same equation as the time-time projection. 

These equations are summarized and contrasted with their electromagnetic counterparts in 
Tabled! Eqs. ([Tj3b)-([Tj6b) are very similar in form to Maxwell Eqs. ([Tj3b)-([Tj6b); they are their 
physical gravitational analogues, since both are the traces and antisymmetric parts of tensors 
{E^tis, Bap} -H- {Eq,/3, Hq,^} that we know, from equations ([Hi) and ^2), that play analogous 
physical roles in the two theories. Note this interesting aspect of the analogy: if one replaces, 
in Eqs. ([U|)- (fT^ . the electromagnetic tidal tensors (Eajs and -Bq/s) by the gravitational ones (Eq,^ 
and Hq,^), and the charges by masses (i.e., density pc and current j° of charge by density p and 
current J°' of mass) one almost obtains Eqs. ([Tj3b)-([Tj6b), apart from a factor of 2 in the source 
term in ([Tj6b) and the difference in the source of Eq. ([9]), signaling that in gravity pressure and 
stresses contribute as sources. This happens because, since E^,^ and Hq,^ are spatial tensors, all 
the contractions with U" present in Eqs. (jlOp and (|12p vanish. That is, in the framework of the 
analogy based on tidal tensors, one starts with Maxwell equations, and by a simple application of 
the analogy one almost ends up with the exact temporal part of Einstein's equations. 

The space-space part of Einstein equations, Eq. ([H?), has no electromagnetic analogue. It 
involves Faf^ and the space-space part of the energy momentum tensor T"^ as a source, none of 
them having an electromagnetic counterpart. This equation also has a fundamental differenc^ with 
respect to the temporal equations ([Tj3b)-([Tj6b) and their electromagnetic analogues: the latter are 
algebraic equations involving only the traces and antisymmetric parts of the tidal tensors, plus the 
source terms; they impose no condition on the symmetric parts. In electromagnetism, this is what 
allows the field to be dynamical, and waves to exist (their tidal tensors are described, in an inertial 
frame, by Eqs. (|128p - ()129p below); were there additional independent algebraic equations for the 
traceless symmetric part of the tidal tensors, and these fields would be fixed. But Eq. (f2T]) . by 
contrast, is an equation for the symmetric parts of the tensors E^^ and Eq,^. It can be split in two 
parts. Taking the trace, and using ([Il3b), one obtains the source equation for ¥a/3- 



M", = 0; (a) F[„^]=0; (b) E[,^] = (c). 



(22) 



Svrp 



(23) 



substituting back in ([H?) we get: 




(24) 



■^We thank Joao Penedones for drawing our attention to this point. 
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This equation tells us that the tensor F"^ is not an extra (comparing with electrodynamics) in- 
dependent object; given the sources and the gravitoelectric tidal tensor E„^, F^^ is completely 
determined by (fM|) . 

A remarkable thing occurs however in vacuum (T"^ = 0, = 0). The Riemann tensor becomes 
the Weyl tensor: Rap-^s = Cap-yS', due to the self duality property of the latter: Ca/s-fS = —*C-kai3'y5, 
it follows Fai3 = — Fq^, and therefore Eq. ([117) identically vanishes. This means that in this case, 
the full Einstein equations and the algebraic Bianchi identity reduce to equations for the traces and 
antisymmetric parts of the the tidal tensors, just like the Maxwell equations. This is summarized 
in Table [3 

Table 2: In vacuum, the full Einstein equations, plus the algebraic Bianchi Identity, reduce to 
equations for the traces and antisymmetric parts of the tidal tensors, just like Maxwell's equations. 

Vacuum 



Maxwell's Equations 




Einstein Eqs. 


+ Bianchi Id. 


^. = 


(Ella) 


F"„ = 


dUlb) 




(I2l2a) 


= 


^2h) 


B^^ = 


da 3a) 


M'^^ = 


(I2l3b) 




(I2l4a) 


= 


(I2l4b) 



These algebraic equations, as well as their general version ([Tl3)-([Tl6) of Table [H have the 
status of constrains for the tidal fields. They are especially suited to compare the tidal dynamics 
(i.e., Eqs. ([Hi) and ([Il2)) of the two interactions, which is discussed in the next section; such 
comparison allows us to learn some fundamental aspects of both interactions, which is exemplified, 
in the context of the dynamics of spinning test particles, in the companion paper But they do 
not tell us about the dynamics of the fields themselves. To obtain dynamical field equations, one 
possible route is to take one step back and express the tidal tensors in terms of gauge fields (such 
as the GEM "vector" fields G, H and the shear K(^a;i3) of the 1+3 formalism of Sec. [Sj the general 
expressions of Einstein equations in terms of these fields is given in Eqs. (7.9) and (7.3) of [19]); but 
it is also possible to write the equations for the dynamics of the tidal tensors (the physical fields); 
that is done not through Einstein equations (jl4p . but through the differential Bianchi identity 
RaT[ni';a] = 0, together with decomposition ([T5|) . and using ([H]) to substitute by the source 
terms. The resulting equations, for the case of vacuum (where {£ai3, ^a/?} = Hq^}), are 

Eqs. (|12ip - (jl22p of Sec.[6]below. One may write as well dynamical equations for the electromagnetic 
tidal tensors, which for the case of vacuum, and an inertial frame, are Eqs. (|128p - (|129p of Sec. 16.11 
however in the electromagnetic case the fundamental physical fields are the vectors E", B"', whose 
covariant field equations are Eqs ([MD, ([58]), ([60|), ([M]) (the tidal field equations ([T28]> - (fT29]l foUow 
trivially from these). 
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2.2 Gravity vs Electromagnetism 

The tensor — whereas Maxwell's equations can be fully expressed in terms of tidal tensors 
and sources, the same is only true, in general, for the temporal part of Einstein's equations. The 
Space-Space part, Eq. (j2ip . involves the tensor F^^, which has no electromagnetic analogue. F^^ 
however is not an additional independent object, as it is completely determined via ()2ip given the 
sources and E^^. In vacuum F^^ = — E^^, and Einstein equations (plus the algebraic Bianchi 
Identity) reduce to equations for the traces and antisymmetric parts of the tidal tensors, analogous 
to their electromagnetic counterparts, cf. Table [2j 

Sources — The source of the gravitational field is the rank two energy momentum tensor T"^, 
the source of the electromagnetic field is the current 4- vector j". Using the projectors ([2|) one 
can split T"^ = pVU^ + 2/i"^ J^'^C/'^) + T^")^'^^, and = pJJ"" + 2h°'^f. Eqs. ([Il3) show that 
the source of Eai3 is pc, and its gravitational analogue, as the source of Eq,^, is 2p + T*^ (p + 2>p 
for a perfect fluid). The magnetic/gravitomagnetic tidal tensors are analogously sourced by the 
charge/mass-energy currents j^''^/J^'^\ as shown by Eqs. ([Tj6). Note that, in stationary (in the 
observer's frame) setups, ^F^/s-'yU'^ vanishes and equations ([Tj6a) and ([Tj6b) match up to a factor 
of 2, identifying j'^'*^ -H- J^^K Eq. (I23p shows that p is the source of Fq,^. Eq. ([Il7), sourced by the 
space-space part T^°'^^^\ as well as the contribution T'^ for ([HSb), manifest the well known fact 
that in gravity, by contrast with electromagnetism, pressure and stresses act as sources of the field. 

Symmetries and time projections of Tidal Tensors — The gravitational and electromagnetic 
tidal tensors do not generically exhibit the same symmetries; also the former tidal tensors are 
spatial, whereas the latter have a time projection (with respect to the measuring them), signaling 
fundamental differences between the two interactions. In the general case of fields that are time 
dependent in the observer's rest frame (that is the case of an intrinsically non-stationary field, 
or an observer moving in a stationary field), Eai3 possesses an antisymmetric part, which is the 
covariant derivative of the Maxwell tensor along the observer's worldline; Eq,^, by contrast, is always 
symmetric. As discussed above, Ea/s is a covariant derivative of the electric field as measured in 
the MCRF; and Eq. ()13b ) is a covariant way of writing the Maxwell-Faraday equation V x E = 
—dB/dt. Therefore, the statement encoded in the equation E[q,^] = 0, is that there is no physical, 
gravitational analogue to Faraday's law of induction (in the language of GEM vector fields of 
Sec. [3l we can say the the curl of the gravitoelectric field G does not manifest itself in the tidal 
forces, unlike its electromagnetic counterpart. See Sec. 13.51 for explicit demonstration). To see a 
physical consequence, let in Eq. ([TJla) — the separation vector between a pair of particles with 
the same q/m and the same 4- velocity C/° — be spatial with respect to U"' {6x°'Ua = 0); and note 
that we can write the spatially projected antisymmetric part of Efj_u can be written in terms of the 
dual spatial vector a^: i^j^^^^j,^] = e^u^ySCx'^U^ . Then the spatial components ([Ula) can be written 
as (using ^(^)(,.) = + E[(^ij.){u)])- 



D'^^xi^^) _ q 



dT^ 



m 



s ' D'^Sx q " ( ) 

E{{n){v))^x'^ + ^fii/fS'^'^U Sx'^ <^ ^ ^ = — E ■ 5x + 6x X a (25) 



the second equation holding in the frame W = 0, and we used the dyadic notation ^ of e.g. |105j . 
From the form of the second equation we see that qa/m is minus an angular acceleration. Using 
relation ([1]), we see that a^' = -B'^^U^^; and in an inertial frame a = dB/dt = —V x E. In the 
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gravitational case, since E^^ = ^(^^u) = ^(fj.)(u): we have 

That is, given a set of neighboring charged test particles, the electromagnetic field "shears" them 
via -E(^jy) , and induces an accelerated rotation via the laws of electromagnetic induction encoded in 
The gravitational field, by contract, only shearflthe particles , since E^^^j = 0. 

Further physical evidence for the absence of a physical gravitational analogue for Faraday's law 
of induction is given in the companion paper [3]: consider a spinning spherical charged body in an 
electromagnetic field; and choose the MCRF; if the magnetic field is not constant in this frame, 
by virtue of equation V x E = —dB/dt a torque will in general be exerted on the body by the 
induced electric field, changing its angular momentum and kinetic energy of rotation. By contrast, 
no torque is exerted in a spinning "spherical" body (i.e., a particle whose multipole moments in a 
local orthonormal frame match the ones of a spherical body in flat spacetime) placed in an arbitrary 
gravitational field; its angular momentum and kinetic energy of rotation are constant. 

There is also an antisymmetric contribution •*^Fai3--yU'^ to B^j^; in vacuum, Eq. ([Tj6a) is a 
covariant form of V x i? = dE/dt; hence the fact that, in vacuum, ^[ap] = 0) means that there is 
no gravitational analogue to the antisymmetric part -Bj^^] (i.e., the curl of B) induced by the time 
varying field E. Some physical consequences of this fact are explored in 0]: Eq. ([Tl6a) implies, 
via ([Tj2a), that whenever a magnetic dipole moves in a non-homogeneous field, it measures a non 
vanishing B^^iS] (thus also Ba/3 / 0) , and therefore (except for very special orientations of the dipole 
moment a force will be exerted on it; in the gravitational case, by contrast, the gravitational 
force on a gyroscope is not constrained to be non-vanishing when it moves in a non-homogeneous 
field; it is found that it may actually move along geodesies, as is the case of radial motion in 
Schwarzsclnldspacetimi , or circular geodesies in Kerr-dS spacetime. 

The spatial character of the gravitational tidal tensors, contrasting with their electromagnetic 
counterparts, is another difference in tensorial structure related to the laws of electromagnetic 
induction: as can be seen from Eqs. (jlOp and (|12p . the antisymmetric parts of the E^p and B^/s 
(in vacuum, for the latter) consist of time projections of these tidal tensors. Physically, these time 
projections are manifest for instance in the fact that the electromagnetic force on a magnetic dipole 
has a non- vanishing projection along the particle's 4- velocity U°', which is the rate of work done 
on it by the induced electric field [Hll], and is refiected in a variation of the particle's proper mass. 
The projection, along U°', of the gravitational force ([Il2b), by its turn, vanishes, and the gyroscope 
proper mass is constant. 



''if the two particles, instead of being free (i.e., moving along geodesies), were connected by a "rigid" rod, the 
symmetric part of the electric tidal tensors would also, in general, torque the rod; hence in such system, we would 
have a rotation, even in the gravitational case, see [T7] pp. 154-155. The same for a quasi-rigid extended body; yet, 
even in this case the effect due to the symmetric parts are very different from the ones arising from electromagnetic 
induction: firstly the former do not require the fields to vary along the particle's worldline, they exists even if the 
body is at rest in a stationary field; secondly, they vanish if the body is spherical, which does not happen with the 
torque generated by the induced electric field, see [1] . 

■^This example is particularly interesting in this discussion. In the electromagnetic analogous problem: a magnetic 
dipole in (initially) radial motion in the Coulomb field of a point charge suffers a force; that force, as shown in [3], 
comes entirely from the antisymmetric part of the magnetic tidal tensor: Bap = B^^^] ; it is thus a natural realization 
of the arguments above that H^/? = in the analogous gravitational problem. 
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2.3 The analogy for differential precession 

Eqs. ([Il2) of Table [1] give the tensors and M^i^ a physical interpretation as the tensors which, 
when contracted with a magnetic/gravitomagnetic dipole vector, yield the force exerted on magnetic 
dipoles/gyroscopes. We will now show that these tensors can also be interpreted as tensors of 
"relative", or "differential", precession for these test particles; i.e., tensors that, when contracted with 
a separation vector Sx^ , yield the angular velocity of precession of a spinning particle (possessing 
magnetic moment) at given point V2 relative to a system of axes anchored to spinning particles, 
with the same ^-velocity, at the infinitesimally close point Vi. Somewhat analogous to the electric 
tidal tensors E^jB and E^^, which, when contracted with 5x^ , yield the relative acceleration of two 
infinitesimally close test particles with the same 4-velocity. 

For clarity we will treat the gravitational and electromagnetic interactions separately. We will 
start by the gravitational problem. Consider a system of Fermi coordinates (see e.g. |100lll01j ) with 
origin along a worldline passing trough the location Pi, and denote by its basis vectors. Assuming 
Mathisson-Pirani spin condition S'^^Uj3 = 0, the spin vector of a gyroscope undergoes Fermi- Walker 
transport, see e.g. [4]. Hence, a gyroscope at Vi, momentarily at rest in this frame (i.e., U = eo) 
does not precess relative to it, as by definition the axes are Fermi- Walker transported therein. 
However at V2, a gyroscope (with the same 4-velocity U") will be be seen to precess relative to 
this frame; the angular velocity of precession is minus the angular velocity of rotation of the basis 
vectors relative to Fermi- Walker transport at V2- Let X" be a vector fixed with respect to the 
basis vectors at V2', it has the covariant derivative: 

^ = — + r^^{V2)u^x^ = rij{V2)x^ = R'j^oX^Sx'^ 

where we used the expression for the connection TQj{V2) = R'' ji^^dx^ , see e.g. |101j . From Eq. (jl]), 
we note that 

which, in the Fermi frame at Vi (orthonormal, where C/* = 0) reads: Rij-yr = —(-ijk * R 
We thus have Tqj{V2) = -e^^f^i^R^^ 5x^ = -e' j,M\6x'- and therefore 



fco 



7r 



= X X 5J1g; = tt^V^' (27) 

Hence, a gyroscope at V2 precesses with angular velocity 6Qg relative to a system of axes fixed to 
guiding gyroscopes at Vi (provided that the gyroscopes have the same 4-velocity V^). This result 
was first obtained in a recent work |102j through a different procedure; we believe the derivation 
above is more clear, and shows that one of the assumptions made in |102| to obtain 5Qg — that 
the gyroscopes at Vi and V2 have the same acceleration — is not necessary, they only need to be 
momentarily comoving (i.e., have the same 4-velocity) in order for (|27p to hold. The differential 
dragging effect, in terms of the eigenvalues of and its associated field lines, and its visualization 
in different spacetimes, are discussed in detail therein. 

Let us turn now to the analogous electromagnetic problem. Consider, in fiat spacetime, a triad 
of orthogonal magnetic dipoles with 4-velocity 17°' at Vi; let (edip)a be an orthonormal tetrad at 
Vi such that (edip)Q = U and the spatial triad (edip)- are axes attached to the dipoles. This will 
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be the reference system. If Mathisson-Pirani condition holds, the spin vector of a magnetic dipole 
(magnetic moment /x" = aS°^, a = gyromagnetic ratio) evolves as (e.g. [1]): 



whose space part reads: 



n c-M 

= S.a-U^' + ae^^p^U'^S'^B^iVi) (28) 

— = 5 X ^{Vi); n{Vi) = aB{Vi) 

i.e., the magnetic dipoles (and spatial triad (e^ip)- attached to them) precess relative to Fermi- 
Walker transported axes with an angular velocity —Vl{Vi). Consider now a test magnetic dipole at 
7^2) momentarily moving with the same 4- velocity C/". It will likewise precesses with an angular 
velocity —0(^2) = —crB{V2) relative to Fermi- Walker transported axes at ■ Now consider the 
global inertial frame (since we are in flat spacetime) momentarily comoving with the two sets of 
dipoles; the Fermi- Walker transport law for a spatial vector X" [U^Xa = 0), is DX /dr = 0, which 
in the inertial frame reads dX/dr = 0. That is, the spatial triads of Fermi- Walker transported 
tetrads do not rotate relative to a momentarily comoving inertial frame. And that effectively means 
that, in flat spacetime, the spatial axes of two Fermi- Walker transported tetrads along two distinct 
curves do not rotate one with respect to another provided that they are momentarily comoving 
(regardless of the acceleration of each of them). Hence, since the points V2 and Vi are taken to be 
momentarily at rest, the magnetic dipole at position V2 precesses with respect to the systems of 
Fermi Walker transported axes at 7^2 or Vi with the same frequency —Q.(J^2), and therefore with 
frequency (511em = ^{'Pi) — ^{^2) relative to the triad (e^ip)- attached to the reference dipoles at 
Vi. Again, since Uf = U^ = U",it follows that 

B^{V2) = B'^iVi) + B^Vi)5x^U^ + 0{5x^) 

where we used a Taylor expansion of S° = F°C/^ around Vi, and noted that ★F" ^^fe^f/^ = 
BZ^6x"' . We obtain thus the relative precession rate 

5n\^u = ^'CPi) - ^\'P2) = -aB\5x^ (29) 

manifesting once more the physical analogy Hq,^ -h- B^p, present also in the forces exerted on 
spinning particles, Eqs. ([Il2), and in the field equations ([Il6), ([Il5) of Tabled! 



3 Gravito- electromagnetic analogy based on inertial fields from 
the 1+3 splitting of spacetime 

This approach has a different philosophy from the tidal tensor analogy of Sec. [21 Therein we 
aimed to compare physical, covariant forces of both theories; which was accomplished through 
the tidal forces. Herein the analogy drawn is between the electromagnetic fields i?", i?" and 
spatial inertial fields G", (i-e., fields of inertial forces, or "acceleration" fields), usually dubbed 
"gravitoelectromagnetic" (GEM) fields, that mimic and in gravitational dynamics. Inertial 
forces are fictitious forces, attached to a specific reference frame, and in this sense one can regard 
this analogy as a parallelism between physical forces from one theory, and reference frame effects 
from the other. 
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The GEM 3- vector fields are best known in the context of linearized theory for stationary space- 
times, e.g. [23\ I24j. where they are (somewhat naively) formulated as derivatives of the temporal 
parts of the linearized metric tensor (the GEM potentials, in analogy with the EM potentials). 
More general approaches are possible, if one observes that these are fields associated not to the 
local properties of a particular spacetime, but, as stated above, to the kinematical quantities of 
the reference frame. In particular, the GEM fields of the usual linearized approaches are but the 
acceleration and vorticity of the congruence of "static" observers u° ~ 5q in the chosen background. 
Taking this perspective, the GEM fields may actually be cast in an exact form applying to arbitrary 
reference frames in arbitrary fields, through a general 1+3 splitting of spacetime. In this section 
we present such an exact and fully general formulation. The spacetime is split with respect to a 
congruence of observers of 4- velocity u"', in a time direction parallel to u'^, and the 3-dimensional 
hyperplanes orthogonal to it, which has the interpretation of the rest space of the observers. Asso- 
ciating to each observer a tetrad e^, such that eg = u, and the spatial triad spans its local rest 
space, in the reference frame thereby defined the 4-dimensional physics is mapped to a (general- 
ized, non-holonomic) curved 3-space where there "live" four spatial fields: G", O", and K^^f^^, 
encoding the temporal part of the connections (rggirg-.,r^.g,r?-.). is the "gravitoelectric field" 
(minus acceleration of the observers), 0" and are, respectively, the rotation of the local tetrads 
relative to local Fermi- Walker transport and the vorticity of the congruence, which together form 
the gravitomagnetic field = 0" + uj". The rank 2 spatial tensor K(^a(3) has no electromagnetic 
analogue, and consists of the shear and expansion of the observer congruence. The treatment herein 
is largely equivalent to the approach in |19l [20] , only it is formulated herein in terms of an arbitrary 
reference frame, not necessarily the "congruence adapted" one (defined by setting = w", see 
below). The main difference (apart from the differences in the formalism) is that in this more 
general formulation the gravitomagnetic field H°' explicitly splits in the two independent parts cj" 
and $7" ; this degree of generality is of use both in the companion paper [3] , and in [63] . 



3.1 The reference frame 

To an arbitrary observer of 4- velocity u" at a given point x", one naturally associates an adapted 
frame (e.g. |19]). which is a tetrad whose time axis is the observer's 4- velocity eg = u and 
whose spatial triad e^ spans the local rest space of the observer. The latter is for now undefined up 
to an arbitrary rotation. The evolution of the tetrad along the observer's worldline is generically 
described by the equation: 

where 0.°^^ is the (anti-symmetric) infinitesimal generator of Lorentz transformations, whose spatial 
part Q'-j = describes the arbitrary angular velocity Q. of rotation of the spatial triad e^ 

relative to a Fermi- Walker transported triad. Alternatively, from the definition of the connection 
coefficients, Ve^e^ = P^.e^, we can think of the components of Vt^^ as some of these coefficients: 

a' ; (30) 

(31) 



" 


00 


= r?- 

Oi 
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Unlike the situation in flat spacetime (and Lorentz coordinates), where one can take the tetrad 
adapted to a given observer as a global frame, in the general case such tetrad is a valid frame only 
locally, in an infinitesimal neighborhood of the observer. In order to define a reference frame over 
an extended region of spacetime, one needs a congruence of observers, that is, one needs to extend 
ti" to a field of unit timelike vectors. A connecting vector X" between two neighboring observers 
in the congruence satisfies 

[u, X] = 4^ n^V^X" - X^V/3n" = 0. (32) 

The evolution of the connecting vector along the worldline of an observer in the congruence is then 
given by the linear equation 

VuX" = (v'^u") Xp. (33) 



The component of the connecting vector orthogonal to the congruence, 

y« = X" + (upX^^ u'^, (34) 

satisfies 

V^y" = K'^'^Yfs + (afsYl^^ u", (35) 

where K"'^ denotes the spatially projected covariant derivative of n"the tensor 

i^"^ = {h'')\{h''f^u^'^ = V^u" + a"?x^ (36) 

The decomposition of this tensor into its trace, symmetric trace-free and anti-symmetric parts 
yields the expansion 9 = Vau", the shear 

= K^^^p) - '^Bqo.p - ^OUaUp (37) 

and the vorticity 

= iSTjo,^] (38) 

of the congruence. It is useful to introduce the vorticity vector 

^ l^a/375^^.^^<5 = _ ^ g«/37^^^^^5. (39) 

cj", according to definition above, yields half the curl of u"; this is in agreement with the convention 
in e.g. [141 119] . but differs by a minus sign from the definition in e.g. [321 186j . Note however that 
for the vorticity tensor w^/? we are using the more general definition given in [32^ [86] , differing from 
a minus sign from the one in [19] (consequently, given by Eq. ([39l) is minus the dual of oJap)- 
The nonvanishing tetrad components of K^p are 

K~.~. = a.. + -e5..+u:... (40) 
These components determine the following connection coefficients: 

K'.^ = V.nj = Vi. = t\.. (41) 
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The remaining temporal connection coefficients (otlier than the ones given in Eqs. (j30p - (j3ip . (j4ip 
above) are trivially zero: 

= -^0 • Ve.eg = -^Vejeg • eg) = 0. 

Each observer of the congruence carries its own adapted tetrad, and to define the reference 
frame one must provide the law of evolution for the spatial triads orthogonal to u°'. A natural 
choice would be Fermi- Walker transport, $7 = (the tetrad does not rotate relative to local guiding 
gyroscopes); another natural choice, of great usefulness in this framework, is to lock the rotation of 
the tetrads to the vorticity of the congruence, Q = u). This has been dubbed in the literature |191l20j 
"co-rotating Fermi- Walker transport" (and the frame thereby defined is "the frame adapted to the 
congruence", which is regarded as the most natural generalization of the non-relativistic concept 
of reference frame, see I37|). This choice is more intuitive in the special case of a shear- free 
congruence, where, as we will show next, the axes of the frame thereby defined point towards fixed 
neighboring observers. Indeed, if X" is a connecting vector between two neighboring observers of 
the congruence and is its component orthogonal to the congruence, we have 

= Y'' + n'-.xK (42) 



Since from (1351) we have 



we conclude that 



V^Y' = K\Y^, (43) 
K = ( ar, + -95 + uj.. - n.A Y^- (44) 



This tells us that for a shear- free congruence (crj-- = 0), if we lock the rotation of the tetrad 
to the vorticity uj of the congruence, VL-fj — '^Ip the connecting vector's direction is fixed on the 
tetrad (and if in addition = 0, i.e., a rigid congruence, the connecting vectors have constant 
components on the tetrad). A familiar example is the rigidly rotating frame in fiat spacetime; in 
the non-relativistic limit, the vorticity of the congruence formed by the rigidly rotating observers 
is constant, and equals the angular velocity of the frame; in this case, by choosing = tj, one is 
demanding that the spatial triads carried by the observers co-rotate with the angular velocity of 
the congruence; hence it is clear that the axes always point to the same neighboring observers. 
For relativistic rotation, the vorticity u is not constant and no longer equals the angular velocity of 
the rotating observers; but the condition = w still ensures that the tetrads are rigidly anchored to 
the observer congruence. Another example is the family of the so-called "static" observers in Kerr 
spacetime, which is very important in this context, because it is this construction which allows one 
to determine the rotation of the frame of the "distant stars" with respect to a local gyroscope, as 
explained in Sec. 13.31 

3.2 Geodesies — "gravitoelectromagnetic fields" 

The equation of geodesies for a test particle of 4-velocity C/", Vu^^" = DU°' /dr = 0, reads, in the 
frame e^: 

w + ([75)2 + (r|. + pL) u^w + rijdjf'u^ = o 
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Substituting pOj) . (j3T]) and dHJ, the spatial part becomes 



dr 



U^G + U X H- a^.uh'. - -eU 
3 ' 3 



(45) 



where DW/dr = W + T\~U^U^ denotes the tetrad components of the spatially projected covariant 

derivative, DU"'/dT = {h^)°'i^DU^ /dr , G = —a is the "gravitoelectric field", and H = Co + 0. is the 
"gravitomagnetic field". These designations are due to the analogy with the roles that the electric 
and magnetic fields play in the electromagnetic Lorentz force, which in this notation reads 

^ = ^(u'^E + Uxb), (46) 
dr rriQ \ J 

with E = E{u) and B = B{u) denoting the electric and magnetic fields as measured by the observers 
w". It is useful to write the GEM fields in a manifestly covariant from: 

(G")° = -Vutx" = -^xV^; (//")" = + 0", (47) 

The gravitomagnetic field (iJ")" consists thus of two parts of diff'erent origins: the angular velocity 
of rotation of the tetrads relative to Fermi- Walker transport (i.e., the local guiding gyroscopes), 
plus the vorticity cj" of the congruence of observers u". If we lock the rotation of the tetrad to the 
vorticity of the congruence: fi" = w", the gravitomagnetic field becomes simply twice the vorticity: 

The last two terms of ()45p have no electromagnetic counterpart; they correspond to the time 
derivative of the spatial metric. In fact, they involve only K(^a(3)^ which is sometimes called the 
second fundamental form of the distribution of hyperplanes orthogonal to u. If this distribution is 
integrable (that is, if there is no vorticity) then is just the extrinsic curvature of the time 

slices orthogonal to u. 

The quotient of the spacetime by the congruence is the "space of observers". In the general case, 
where 7^ 0, there is no natural metric one can associate to this space, as the distance between 

observers, unlike the quotient, depends on the time slices. In the special case that the congruence 
is rigid, one naturally associates to this space the metric measuring the constant distance between 
neighboring observers, which is is obtained by lifting tangent vectors on the quotient to tangent 
vectors orthogonal to the corresponding curve. On the other hand, if we have K[ai3) 7^ 0, but 
o;"^ = 0, the congruence is hypersurface orthogonal, and it is natural to define a (time-dependent) 
3-D metric on those hypersurfaces. 

It is worth noting that Eq. (|45p is the most general description of the geodesic motion possible in 
terms of electromagnetic- like fields (and kinematical quantities with no electromagnetic analogue) ; 
the corresponding results presented in the more popular linearized theory [241 123] or Post-Newtonian 
[691 [mils] approaches are special cases of this equation (e.g. linearizing Eq. ([15]) above, one obtains 
Eq. (2.5) of [25]; further specializing to stationary fields, one obtains e.g. (6.1.26) of [23]). 

3.2.1 Stationary fields — "Quasi-Maxv^ell" formalism 

If one considers a stationary spacetime, and a frame where it is explicitly time- independent (i.e., 
a congruence of observers tangent to a time-like Killing vector field, which necessarily means 
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that the congruence is rigid |118j ) , the last two terms of Eq. (j45p vanish and the geodesic equation 
becomes formally similar to the Lorentz force (06]): 



DU ud (ijOg + U X H) (48) 



dT 

The line element of a stationary spacetime is generically described by: 

ds^ = -e^'^idt - Aidx'f + -Hjdx'dx^ (49) 

with A, 'Jij time-independent, jij is a spatial metric, not flat in general, that measures the 
constant distance between stationary observers, as measured by the Einstein light signaling pro- 
cedure |77] . This is the metric one can naturally associate with the quotient space in the case a 
rigid congruence, as discussed above. The GEM fields measured by the static observers (i.e. the 
observers of zero 3- velocity in the coordinate system of I49p are related with the metric potentials 
by dU: 

G = -V$; H = e*V x A, (50) 

with V denoting the covariant differentiation operator with respect to the spatial metric jij. The 
formulation (|50p of GEM fields applying to stationary spacetimes is the most usual one; it was 
introduced in [77], and further worked out in e.g. [El [HI (TT] (93] [92], and is sometimes called the 
"Quasi-Maxwell formalism". 

3.3 Gyroscope precession 

One of the main results of this approach is that within this formalism, the equation describing the 
evolution of the spin- vector of a gyroscope in a gravitational field (i.e., the Fermi- Walker transport 
law): 

— ^ = 5,a^ [7° , (51) 
dr 

takes, when expressed in a local orthonormal tetrad comoving with the test particle, a form exactly 
analogous to precession of a magnetic dipole under the action of a magnetic field. This analogy is 
more general than the one for the geodesies described above, it holds for arbitrary fields, and the 
observer does not have to be stationary (i.e. its worldline does not have to be tangent to a Killing 
vector); the only condition is to be comoving with the gyroscope. 

Let C/° be the 4- velocity of the gyroscope; in a comoving orthonormal tetrad Ca, = (5?, and 
also = 0; therefore, Eq. (|5ip reduces in such frame to: 



dr dr 

This result is somewhat obvious; note that we are just saying that, relative to this frame, gyroscopes, 
which are objects that "oppose" to changes in direction, and determine the the local "compass of 
inertia", are seen to "precesses" with an angular velocity that is simply minus the angular velocity 
of rotation of the frame relative to Fermi- Walker transport. Now, for a congruence adapted frame, 
f1 = w, this becomes: 

%^\S.H (52) 
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Table 3: The gravito-electromagnetic analogy based on inertial GEM fields. 



Electromagnctism 



Gravity 



Lorentz Force: 



Geodesic Equation {}j = q . 



cLt mo V / 



(Ella) ^ = C/° 
dr 



U°G + U xH- a'-.U^e-- - -0U 
3 * 3 



Plb) 



Precession of magnetic dipole: 








Gyroscope "precession": 




DS ^ ^ 
dr ^ 


(ISl2a) 






dr 


(I3l2b) 


Stationary fields, rigid, congruence 


adapted fr 


ame: D, 


= u) = H /2 (Quasi-Maxwell formalism) 


Force on magnetic dipole; 








Force on gyroscope: 




^ - V(B • a/) - ^^(V • B) - 1(^1 • H)E 


il3a) 


DP 


1 
2 


V(i? • S) - S{\/ ■ H) - 2{S-H)G 


(I33b) 


Maxwell Source Equations 








Einstein Equations 














• Time Component! 






• Time-Time Component: 




V -E^ i-Kpc + H ■ B 


(I3l4a) 


V 


■G 


= -4^(2p + r°) + G2 + iF2 


(1314b) 


• Space Components: 






• 


Time-Space Components: 




\/xB^GxB + Attj 


(01 5a) 




VxH = 2GxH- 16nJ 


(I35b) 








• 


Space-Space Component: 




No electromagnetic analogue 




y^G, - 




-t- 71116) 


Bianchi Identity 






Algebraic Bianchi Identity 














• Time Component: 




• Time-Time (or Space-Space) Component: 




V ■ B = -H ■ E 


dSlTa) 






V ■ H = -H ■ G 


(Il7b) 


• Space Components: 






• 


Space-Time Components: 




y X E = G X E 


ilSa) 






V X G = 


(I38b) 



Thus, the "precession" of a gyroscope is given, in terms of the gravitomagnetic field H, by an 
expression identical (up to a factor of 2) to the precession of a magnetic dipole under the action of 
a magnetic field B, cf. Eq. (p8l) : 

D5/fiT = flxB. (53) 

This holds for arbitrary fields, and hence the result obtained for of weak fields in [25] (that the 
analogy holds even if the fields are time dependent) is just a special case of this principle. However 
important differences should be noted: whereas in the electromagnetic case it is the same field B 
that is at the origin of both the magnetic force q{U x B) in Eq. (|46|) and the torque fl x B on the 
magnetic dipole, in the case of the gravitomagnetic force U x H it has, in the general formulation, a 
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different origin from gyroscope "precession", since the former arises not only from the rotation of 
the frame relative to a local Fermi- Walker transported tetrad, but also from the vorticity a; of the 
congruence. Even in the case = uj there will be a relative factor of 2. In this sense, one can say 
that the Lense-Thirring effect detected in the LAGEOS satellite data |41j (and which will also be 
subject of experimental detection by the upcoming LARES mission [l3]), measuring H from test 
particle's deflection, is of a different mathematical origin from the one which was under scrutiny 
by the Gravity Probe B mission [42], measuring il. from gyroscope precession, the two being made 
to match by measuring both effects relative to the "frame of the distant stars" (see below) . 

Another obvious difference is the presence of a covariant derivative in (|53p , and a simple deriva- 
tive in (I52p . signaling that i? is a physical field, and H a mere artifact of the reference frame (which 
can be anything, depending on the congruence of observers one chooses), that can be made to 
vanish by choosing a vorticity-free congruence. Note however that this does not mean that H is 
necessarily meaningless; indeed it has no local physical significance, but it can tell us about frame 
dragging, which is a non-local physical effect. That is the case when one chooses the so-called 
"frame of the distant stars", a notion that applies to asymptotically flat spacetimes. In stationary 
spacetimes, such frame is setup as follows: consider a rigid congruence of stationary observers such 
that at infinity it coincides with the asymptotic inertial rest frame of the source — the axes of the 
latter define the directions fixed relative to the distant stars. If the spacetime is also axisymmet- 
ric (such as the Kerr spacetime), this congruence asymptotically coincides with the zero angular 
momentum observers (ZAMOS, |117l [28| [T7]). These observers are interpreted as being "at rest" 
with respect to the distant stars (and also at rest with respect to the asymptotic inertial frame 
of the source); since the congruence is rigid, it may be thought as a grid of points rigidly fixed 
to them. For this reason we dub them "static observers'[§. This congruence fixes the time axis 
of the local tetrads of the frame. Now if we demand the rotation of the local tetrads (relative 
to Fermi- Walker transport) to equal the vorticity uj of the congruence, Eq. (j44p shows that the 
connecting vectors between different observers are constant in the tetrad; in other words, each local 
spatial triad e~- is locked to this grid, and therefore has directions fixed to the distant stars. Hence, 
despite having no local meaning, the gravitomagnetic field H = 2Q = 2io describes in this case a 
consequence of the frame dragging effect: the fact that a torque free gyroscope at finite distance 
from a rotating source precesses with respect to an inertial frame at infinity. This is a physical 
effect, that clearly distinguishes for instance the Kerr from the Schwarzschild spacetimes, but is 
non-local (i.e., it cannot be detected in any local measurement; only by locking to the distant stars 
by means of telescopes). It should be noted however that the relative precession of two neighboring 
(comoving) system of gyroscopes is locally measurable and encoded in the curvature tensor (more 
precisely, in the gravitomagnetic tidal tensor M^fs, as discussed in Sec. 12.31 

®In the case of Kerr spacetime, these are the observers whose worldhnes are tangent to the temporal KiUing vector 
field ^ = d/dt, i.e., the observers of zero 3- velocity in Boyer Lindquist coordinates. This agrees with the convention 
in [1171128) . We note however that the denomination "static observers" is employed in some literature (e.g. [10311104] ) 
with a different meaning, where it designates hypersurface orthogonal time-like Killing vector fields (which are rigid, 
vorticity-free congruences, existing only in static spacetimes). 



24 



3.4 Field equations 

The Einstein field equations and the algebraic Bianchi Identity, Eqs. (|14p . can be generically written 
in this exact GEM formalism — i.e., in terms of G, 0,, Co and Kf^^i^)- These equations are compared 
with the analogous electromagnetic situation: Maxwell's equations in an arbitrarily accelerated, 
rotating and shearing frame. The latter will be of use also in Sec. [6l We will also consider the 
special case of stationary spacetimes (and rigid, congruence adapted frames therein, Q, = Co = H/2), 
where we recover the Quasi-Maxwell formalism of e.g. |14[ [Ml [771 [TBI 1231 In this case there 
is a closer similarity with the electromagnetic analogue — Maxwell's equations in an arbitrarily 
accelerated, rotating rigid frame (acceleration a = —G, vorticity Co = Q = H/2, cf. Sec. 13. ip . 

3.4.1 Maxwell equations for the electromagnetic fields measured by an arbitrary con- 
gruence of observers 

Using decomposition ([T|), we write Maxwell's Eqs. ([6]) in terms of the electric and magnetic fields 
(£'")" = F'^pU^ and (i?")" = -kF"-pU^ measured by the congruence of observers of 4- velocity u". All 
the fields below are measured with respect to this congruence, so we may drop the superscripts: 
(£"")" = E'", (i?")" = i?". The time projection of Eq. ^) with respect to (see point [2] of Sec. 
II. ip reads: 

E^.p = Anpc + E"a„ + 2w„S" (54) 

To write this equation in the tetrad components of the reference frame given in Sec. 13.11 we note 
that, for a given vector A": 

= + + A'rI. =V-A + A.a, (55) 

where V denotes the spatially projected covariant derivative of a tensor A^''-'"^": 

VaA^'-^- = {K')\{K')%,..\h''f;VxA''^---P-. (56) 

Mathematically, this covariant derivative is a connection in the vector bundle determined by the 
distribution of hyperplanes orthogonal to the congruence. For the special case of a rigid congruence, 
discussed in Sec. 13.2.11 it becomes a covariant derivative with respect to the natural spatial metric 
defined in the quotient space. Therefore, we have 

V ■E = Airpc + 2C0-B, 
or, in the congruence adapted frame {Co = Q = H/2), 

V ■E = iirpc + H-B . (57) 
Analogously, for the time projection of ([lb) we get 

bI^.^ = B''aa-2Lof'E^ (58) 

which in the tetrad becomes 

V ■ B = -2Cj ■ E, 
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or 

V ■B = -H-E (59) 
in the congruence adapted frame. The space projection of Eq. ([6^) reads: 

e^^^^Bp,^ = ^^^^ - K^^I^^E^ + - e^^^w^^T + e^^^S'^a^ + 47rj<"> (60) 

where the index notation (fi) stands for the spatiaUy projected part of a vector: V^^^ = h^Vy. Here 
we used e^^" = e^^'^'^Ua, and the fact that the spatially projected covariant derivative of a spatial 
vector equals its Fermi- Walker derivative: 

DA^ _ DfA- 

T = Tu denotes herein the proper time along a worldline of tangent vector The tetrad compo- 
nents of (j60p in the frame defined in Sec. 13.11 read: 

(V X BY = K^^^'^E'. + 9E' - (uj x EY + (G x BY + Airf. (61) 

ar ^ 

Using the connections ()4ip . we note that the covariant derivative of a vector A"' can be written, in 
the tetrad, as: 

choosing the congruence adapted frame = wj~), and substituting in ([6T]) . we obtain: 

VxB = — + GxB + ATrj- K^'^^E-.e- + BE; (63) 

in a rigid frame [K^^^^ = ^ = 0) where the fields are time-independent {dE/dr = 0), this yields 
Eq. da 5a) of Tabled 

The space projection of dlb) is 

e'^^^Ep,^ = -^^^ + K^^f^^Bp - OB'^ + e^w^S^ + e'^'^'E^a,, (64) 
which analogously becomes, in the congruence adapted frame, 

V xE = — - + GxE + K^'^^B'.e. - OE. (65) 
ar 

If the frame is rigid, and the fields time-independent, this yields Eq. ([3j8a) of Table [3l 
3.4.2 Einstein equations 

We start by computing the tetrad components of the Riemann tensor in the frame of Sec. 13.11 

%6i = -^^^ + G^^-^K., + K^. + dji.,-K'ji., (66) 

%f = + + (68) 
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In expressions above we kept independent of uj, so that they apply to an arbitrary orthonormal 
tetrad field. Here 



R, 



p/ _|_ p£ 



(69) 



ij k jk,i ik,j im jk jrh ik ij rhk 

denotes the Riemann tensor of the distribution of hyperplanes orthogonal to the congruence, that 
is, the curvature of the connection V on the vector bundle determined by the distribution. is 
anti-symmetric in the first and second pair of indices, but in general does not satisfy the algebraic 
Bianchi identities. 

We shall now compute the tetrad components of the Ricci tensor, but now specializing to 



congruence adapted frames: 



eijkH^/2, so that the Ricci tensor comes in 



terms of the three GEM fields: G, H and -f^(ij)- These read 



00 



R: 



O'i) 



+ -{Vy,H\-{GxH). 



1 

+ 2 



t J {ii) J 



(ij) 



(70) 
(71) 

(72) 



where Hj-j = enkH^ is the dual of H, and R— = R^~-~ is the Ricci tensor of the distribution of 
hyperplanes orthogonal to the congruence; this tensor is not symmetric in the general case of a 
congruence possessing both vorticity and shear. Using T^^ = p and T*^* = J*, the time-time, time- 
space, and space-space components of the Einstein field equations with sources, Eq. (fllb ). read, 
respectively: 



Svr T; 



V-G 

V X H 



2 V " 



-47r(2p + r"„) + G^ + ^H^ 
-167rJ + 2G X H + 2W - 2V',K^''>e. 



dr 



+- 



dr ^ ^ ^ « J W J ^ (b) 



(73) 
(74) 

(75) 



Eqs. (I73]l - (j74]) are the gravitational analogues of the electromagnetic equations (i57|) and (j63ll . re- 
spectively. Eq. ()75p has no electromagnetic counterpart. 

As for the the algebraic Bianchi identities (fT^ ). using the time-time (equal to space- 

space), time-space and space-time components become, respectively: 



V-H 

V X G 



-G-H 

dr J * 



(76) 

(77) 
(78) 



Eqs. (j76p -(|77p are the gravitational analogues of the time and space projections of the electromag- 
netic Bianchi identities, Eqs. ([59|) - ([65]) . respectively. Eq. ([78]) has no electromagnetic analogue. 
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This equations states that, when the observer congruence has both shear (including expansion) and 
vorticity, Rijki does not obey the algebraic Bianchi identities for a 3D curvature tensor. Note this 
remarkable aspect: all the terms in the Maxwell equations (|57p . ()59p and (|65|) have a gravitational 
counterpart in ([73]) . ([76l) and ([77j) . respectively, substituting {E,B} — )• {G,H} and up to some 
numerical factors. As for (j63p . there are clear gravitational analogues in (j74p to the terms G x B 
and the current 47rj, but not to the remaining terms. It should nevertheless be noted that, as 
shown in Sec. [5] below, in the Post-Newtonian regime (or in the "GEM limit" of linearized the- 
ory), the term 2V6 of (j74p embodies a contribution analogous to the displacement current term 
dE/dr of (j63p . The gravitational equations contain, as one might expect, terms with no parallel in 
electromagnetism. 

Most of the differing terms involve the total shear tensor K^^'-y The similarity thus get closer 
if we take the "Quasi-Maxwell" regime, i.e., stationary fields, and a frame adapted to a rigid con- 
gruence of stationary observers: Kij = K^^j]^ = ujij. In this case the expressions above for the 
Riemann and Ricci tensors become the ones given in [T3]. The gravitational field equations (jl4p . 
i.e., (|73p -(|78p. in this regime, are given in TableEJ Therein we drop the hats in the indices, for the 
following reason: as discussed in Sec. 13.21 in this regime there is a natural 3-D metric 7jj on the 
quotient space (measuring the fixed distance between neighboring observers); we thus interpret the 
spatial fields G and H as vector fields on this 3-D Riemannian manifold. The spatially projected 
covariant derivative operator V becomes the covariant derivative of 7ij (as T*-^ =(^^ F*^, i.e. the 
4-D spatial connections equal the connection coefficients for 7jj), and Rij its Ricci tensor, which is 
symmetric (contrary to the general case). The equations in this "Quasi-Maxwell" regime exhibit a 
striking similarity with their electromagnetic counterparts, Eqs. (I3l4a)-([3j8a) of Tabled in spite of 
some natural differences that remain — numerical factors, the source and terms in ([3j4b) with no 
electromagnetic counterpart. We note in particular that, by simply replacing {i?, i?} — t- {G,H} in 
([3l5a)-(l3l8a), one obtains, up to some numerical factors, Eqs. ([3l5b), ([3l7b)-([3l8b). Of course, the 
electromagnetic terms involving products of GEM fields with EM fields, are mimicked in gravity 
by second order terms in the gravitational field. This is intrinsic to the non-linear nature of the 
gravitational field, and may be thought of as manifesting the fact that the gravitational field sources 
itself. Note in this context that the term 2G x H = —16-kpg in Eq. ([31 5b), sourcing the curl of 
the gravitomagnetic field, resembles the electromagnetic Poynting vector pem = E x B/Att] and 
the contribution G^ + H"^ /2 = —AnpG in Eq. (I3j5a), sourcing the divergence of the gravitoelectric 
field, resembles the electromagnetic energy density pem = {E'^ + B'^)/8tt. For these reasons pc 
and pg are dubbed in e.g. [92l [931 IM] gravitational "energy density" and "energy current density", 
respectively. It is interesting to note that, in the asymptotic limit, pG corresponds to the time- 
space components of the Landau-Lifshitz [77] pseudo-tensor t^'^ [45] . One should however bear in 
mind that, by contrast with the electromagnetic counterparts, these quantities are artifacts of the 
reference frame, with no physical significance — at least from a local point of view (see related 
discussion in Sec. 16. ip . 

Finally, it is also interesting to compare Eqs. (|73p -(j78 p with the tidal tensor version of the same 
equations, Eqs. (fT9]) - (j22|) . 
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3.5 Relation with tidal tensor formalism 

The analogy based on the gravito-electromagnetic fields G and H is intrinsically different from 
the gravito-electromagnetic analogy based on tidal tensors introduced in [1]; the latter stems from 
tensor equations, whereas the former are fields of inertial forces, i.e., artifacts of the reference 
frame. A relationship between the two formalisms exists nevertheless, as in an arbitrary frame 
one can express the gravitational tidal tensors in terms of the GEM fields, using the expressions 
for the tetrad components of Riemann tensor Eqs. (|66|) -(|67| ) . This relationship is in many ways 
illuminating, as we shall see; it is one of the main results in this work, due the importance of using the 
two formalism together in practical applications, to be presented elsewhere [63]. Herein expressions 
are to be compared with the analogous electromagnetic situation, i.e., the electromagnetic tidal 
tensors computed from the fields as measured in an arbitrarily accelerating, rotating, and shearing 
frame (in fiat or curved spacetime). 

We start by the electromagnetic tidal tensors; from the definitions of E^p and B^jB in Tabled! 
it follows that 

Using decompositions ([H), and Eq. (|62]) . we obtain the tetrad components (E^-- = B^~- = 0): 

= V.E.-eJ^B^Ki. (79) 

By = V.B. + eJ^E^Ky (80) 
dE 

Elo = -^ + i(^xE). + {GxB). (81) 
Bio = ^ + {nxB).-{GxE). (82) 
or, using Kij = ujij + -^^(jj), and choosing a congruence adapted frame {uj = Q = H/2), 



= V'.B'. + - 



B ■ HS'.'. - B .H-. 



(83) 



E ■ H5y - E.H, 



^ ^ , ^l^^E^K,^, (84) 

dE', 1 



Elo = '^+^^{HxE). + {GxBy. (85) 



dr 2 

= ^ + ^(^xS).-(Gx^)5 (86) 

Let us compute the gravitational tidal tensors. From the definitions of Eq,/3 and Hq,^ in Table [U 
and using the tetrad components of the Riemann tensor, Eqs. ()66p -(f67I). we obtain (Eoq, = "^ab ~ 

% = -^^G^ + GlG^-^K.. + K^. + d.K:^-K^.K^ (87) 
% = -V^^J + %V.c^ + 2G.^. + e>V^i^(.^) 
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For a congruence adapted frame these expressions become: 



2 '■''^ dr J"^ 



[VjH, + iG-H)-f,j-2GjHi 



(89) 
(90) 



In (|90p we substituted V ■ H = G ■ H using Eq. (|77p. Note the formal similarities with the 
electromagnetic analogues (i83]) - (i84|) . All the terms present in Eij and Sjj, except for the last 
term of the latter, have a correspondence in their gravitational counterparts Ejj, Hjj, substituting 
{E^B} — )• — {G, ff} and correcting some factors of 2. However, the gravitational tidal tensors 
contain additional terms which (together with the differing numerical factors) encode the crucial 
differences in the tidal dynamics of the two interactions. The third and forth terms in (1830 . which 
have no electromagnetic analogue, have the role of canceling out the antisymmetric part of V-Gj, 
that is, canceling out the contribution of the curl of G to the gravitoelectric tidal tensor, as can be 
seen from Eq. (I77p . In particular, note that the term —dH^/dr, "inducing", via Eq. (I77p . the curl 
of G (which might lead one to predict gravitational induction effects in analogy with Faraday's law 
of electromagnetism) , is being subtracted in (I89p . meaning that the curl of G does not translate in 
the physical, covariant forces. For instance, it does not induce rotation in a set of free neighboring 
particles (see Eq. ([26P above and discussion therein), nor does it torque an extended rigid body, as 
shown in the companion paper [3]. 

There are some interesting special regimes where the relation between the tidal tensor and 
the inertial fields formalism becomes simpler. One is the "Quasi-Maxwell" regime, i.e., stationary 
spacetimes, and a frame adapted to a rigid (i.e., shear and expansion-free) congruence of stationary 
observers. The gravitational tidal tensors as measured in such frame can be expressed entirely in 
terms of the gravito-electric G and gravitomagnetic H fields; the non-vanishing components are: 



E. 



Mi 



-SI jGi + GiGj + ^ 
1 



1 



2 L 



H'ji, - HjHi 
VjHi + (G • H)j^, - 2GjH, 



(91) 
(92) 



The hats in the indices of these expressions are dropped because, since in the Quasi-Maxwell regime 
there is a natural metric ^ij associated to the quotient space (see Sec. l3.2p . we express these tensors 
in terms of an arbitrary (possibly coordinate) basis in the quotient space (as we did in Sec. 13.2.11 
and in Table [3]), instead of tetrad components. 

The non-vanishing components of the electromagnetic tidal tensors are, under the same condi- 
tions. 



Ei. 



Bij 



B-H^^ 



BjHi 



E ■ H-fij - E.Hi 



(a) 



(a) 



E, 



iO 



iO 



]^{HxE)i + {GxB)i (b) (93) 
]^{H ^B),-{GxE)i (b) (94) 



Thus again, even in the stationary regime, the electromagnetic tidal tensors have non-vanishing 
time components, unlike their gravitational counterparts. The spatial parts, however, are very 
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similar in form; note that replacing {E,B} — )■ —{G,H/2} in (j94p . the time components vanish, 
and one almost obtains the space part (jSSp . apart from the factor of 2 in the third term; and 
that a similar substitution in (j93p almost leads to (|9ip . apart from the term GiGj, which has no 
electromagnetic counterpart. The gravitational and electromagnetic tidal tensors are nevertheless 
very different, even in this regime; namely in their symmetries. Eij is not symmetric, whereas Ejj 
is (the second and third terms in (j9ip are obviously symmetric; and that the first one also is can 
be seen from Eq. ([3l8b) of Table [3]). As for the magnetic tidal tensors, note that, by virtue of 
Eq. ([3l5b), the last term of (i92]l ensures that, in vacuum, the antisymmetric part -^[jy] (i.e., the 
curl of H) is subtracted from Hi-j in (j88p . thus keeping Mij symmetric, by contrast with Bij. This 
can be seen explicitly by noting that in vacuum (j92p can be put in the equivalent form: 

Mij = -- Hi-j - H^i-j-^ + (G • H)jij - IG^Hi^ 

where we used = 2Gyffj], as follows from Eq. ([3l5b). 

Another interesting regime to consider is the weak field limit, where the non-linearities of the 
gravitational field are negligible, and compare with electromagnetism in inertial frames. From 
Eqs. (I83p -(|86p. the non- vanishing components for the electromagnetic tidal tensors measured by 
observers at rest in an inertial frame are: 

Eij = Eij] EiQ = Bij = Bij BiQ = , 

i.e., they reduce to ordinary derivatives of the electric and magnetic fields. The linearized gravita- 
tional tidal tensors are, from Eqs. (I89p -()90p: 

1 dH^ d 1 

~ -G^j + ^^ijk^ - ^^feO; (a) ^ij - - 2^^J + ^^'^^{3rn),l ■ (b) (95) 

Thus, even in the linear regime, the gravitational tidal tensors cannot, in general, be regarded as 
derivatives of the gravitoelectromagnetic fields G and H. Noting, from Eq. (jl07p below, that -ft'(ij) 
is the time derivative of the spatial metric, we see that only if the fields are time independent in 
the chosen frame, do we have Ejj ~ —Gij, Mij ~ —^Hij. 

3.6 Force on a gyroscope 

In the framework of the 1+3 formalism, there is also an analogy [H] relating the gravitational force 
on a gyroscope and the electromagnetic force on a magnetic dipole. This is an analogy different 
from the one based on tidal tensors, and not as general. We start with equations ([Il2) of Table [H 
which tell us that the forces are determined by the magnetic/gravitomagnetic tidal tensors as seen 
by the particle. For the spatial part of the forces, only the space components of the tidal tensors, 
as measured in the particle's proper frame, contribute. Comparing Eqs. (j84p and ()90p . which 
yield the tidal tensors in terms of the electromagnetic/gravitoelectromagnetic fields, we see that a 
close formal analogy is possible only when = in the chosen frame. Thus, a close analogy 

between the forces in this formalism can hold only when the particle is at rest with respect to 
a congruence for which Kt^^p) = Oj that is, a rigid congruence. The rigidity requirement can be 
satisfied only in special spacetimes |118| : it is ensured in the "Quasi-Maxwell" regime — that is, 
stationary spacetimes, and congruences tangent to time-like Killing vector fields therein. 
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Let us start by the electromagnetic problem — a magnetic dipole at rest in a rigid, but arbitrarily 
accelerating and rotating frame. Since the dipole is at rest in that frame, /i° = (0,/i*); hence the 
spatial part of the force is F^j^ = B^^Hj. Substituting ([Mb ) in this expression yields the force 
exerted on the dipole, in terms of the electric and magnetic fields as measured in its proper frame: 



Fem = V(S • /2) + ^ \^ll{E ■H)-{fI- H)E 



Using H ■ E = —V • -B, cf. Eq. ([31 7a), we can re-write this expression as 



(96) 



Fem = V{B-fi)-- [/2(V ■B) + {fl- H)E 



(97) 



Consider now the analogous electromagnetic situation: a gyroscope at rest (i.e., with zero 3- 
velocity, f/* = 0) with respect to stationary observers (arbitrarily accelerated and rotating) in a 
stationary gravitational field; from Eqs. ([Il2b) and (|88p . the force exerted on it is given by: 



^ 2 



V{H ■ S) + S{G ■ H) - 2(5 • H)G 



From Eq. ([3l7b) we have G ■ H = —V ■ H; substituting yields 



(98) 



Fg = - 
2 L 



V(i? • 5) - 5(V • H) - 2(5 • H)G 



(99) 



Note that replacing {il,E,B} — )■ {S,G, H /2} in Eq. ([96]) . one almost obtains ([98]) . except for a 
factor of 2 in the last term. The last term of (j98p -(|99p. in this framework, can be interpreted as 
the "weight" of the dipole's energy [H]. It plays, together with Eq. dSlSb), a crucial role in the 
dynamics, as it cancels out the contribution of the curl of H to the force, ensuring that, in the 
tidal tensor form ([Tl2b), it is given by a contraction of 5" with a symmetric tensor Hq^. This is 
discussed in detail in Sec. 13.51 This contrasts with the electromagnetic case, where the curl of B is 
manifest in Bai3 (which has an antisymmetric part) and in the force. The expression (I99p was first 
found in [14J (where it was compared to the force on a magnetic dipole as measured in the inertial 
frame momentarily comoving with it, case in which the last two terms of [971 vanish): herein we add 
expression (|97p , which is its electromagnetic counterpart for analogous conditions (the frame where 
the particle is at rest, which can be an accelerating and rotating), and shows that the analogy is 
even stronger. 



4 "Ultra-stationary" spacetimes 

Ultra-stationary spacetimes are stationary spacetimes admitting rigid geodesic time-like congru- 
ences. In the coordinate system adapted to such congruence, the metric is generically obtained by 
taking $ = in Eq. (jig]) , leading to: 

ds^ = -(^dt- Ai{x'')dx'^ ^ + -fij{x^)dx'dx^ . (100) 

Examples of these spacetimes are the Som-Raychaudhuri [B^, Van-Stockum interior solution [89], 
and Gddel [90] spacetimes; see [2] for their discussion in this context. This is an interesting class 
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of spacetimes in the context of GEM, due to the close similarity with electrodynamics, which was 
explored in an earlier work [1] by one of the authors: 1) they are exactly mapped \68\ [1], via the 
Klein-Gordon equation, in curved 3-spaces with a gravitomagnetic field; 2) their gravitomagnetic 
tidal tensor is linear [1] (just like in the case of electromagnetism) , and, up to a factor, matches 
the covariant derivative of the gravitomagnetic field. A link between these two properties was 
suggested thereir0; however, the non-vanishing gravitoelectric tidal tensor (while no gravitoelectric 
part is present in the map) was a question left unanswered. Herein, putting together the knowledge 
from the tidal tensor with the inertial force formalisms (Sees. [2] and [3|), we revisit these spacetimes 
and shed some new light on these issues. 

Eqs. (I50p yield the GEM fields corresponding to the frame adapted to the rigid geodesic congru- 
ence, u°' = 6q. They tell us that the gravito-electric field vanishes: G = 0, which is consistent with 
the fact that no electric field arises in the mapping above; the gravitomagnetic field H is linear in 
the metric potentials: 

H = V xA. (101) 

These properties can be interpreted as follows. The fact that G = means that the metric is 
written in a frame corresponding to a congruence of freely falling observers (as their acceleration 
a°' = —G°' is zero); the very special property of these spacetimes (unlike the situation in general, 
e.g. the Kerr or Schwarzschild spacetimes) is that such congruence is rigid, i.e. has no shear, allowing 
the metric to be time independent in the coordinates associated to that frame. The gravitomagnetic 
field, on the other hand, does not vanish in this frame, which means in this context (since the frame 
is congruence adapted, see Sec. 13. ip . that the congruence has vorticity. The equation of motion for 
a free particle in this frame, cf. Eq. (j48p . reduces to : 

^ = U^U X H , (102) 
dr 

similar to the equation of motion of a charged particle under the action of a magnetic field; and 
since H is a linear function of the metric, the similarity with the electromagnetic analogue is indeed 
close. 

Let us now examine the tidal effects. This type of spacetimes have a very special property: the 
gravitomagnetic tidal tensor measured by the observers u" = 5q is linear in the fields (and thus 
in the metric potentials), cf. Eq. (|88p . and, just like in the electromagnetic analogue, it is given by 
the covariant derivative of H with respect to ^ij: 

Mij = -^Hi,j = -e^\Ak-ij . (103) 

(Hoj = Hqo = Eljo for these observers) where semi-colons denote covariant derivatives with respect 
to -jij. This reinforces the similarity with electromagnetism. The gravitoelectric tidal tensor is. 



'^In an earlier work by one of the authors [T] [2] (to whom the exact GEM fields analogy of Sec. [3] was not yet 
known), it was suggested that the above mapping could be interpreted as arising from the similarity of magnetic 
tidal forces manifest in relations (|103|l . It seems, however, to be much more related to the analogy based on GEM 
"vector" fields manifest in Eqs. (|101|) and ([102]). Even though the exact correspondence (|103|l reinforces in some 
sense the analogy, tidal forces do not seem to be the underlying principle behind the mapping, since: i) in it there 
is no electromagnetic counterpart to the non- vanishing gravitoelectric tidal tensor E^/?; ii) the Klein-Gordon Eq. 
□ $ = m'^$ and the Hamiltonian in Sec. IV of [T] is for a (free) monopole particle, no tidal forces contribute. It is 
thus expected to reveal coordinate artifacts such as the fields G, H, and not physical tidal forces. 
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however, non-zero, as seen from Eq. ([87]): 

Kij = ^(^H\j-HjHi^ , (104) 

even though G = 0; and that should not be surprising, for the following reasons: i) it is always 
possible to make G vanish by choosing freely falling observers (this is true in an arbitrary spacetime), 
but that does not eliminate the tidal effects, as they arise from the curvature tensor; ii) in the present 
case of ultrastationary spacetimes, Ec,/? is actually a non-linear tensor in H; which merely reflects 
the fact that, except on very special circumstances, E^^ (unlike its electromagnetic counterpart) 
cannot be thought as simply a covariant derivative of some vector field G. 

The tidal tensor (|104|) exhibits other interesting properties. It vanishes along the direction of 
the gravitomagnetic field H": if X" is a vector spatial with respect to n° (i.e., X"ua = 0) which is 
parallel to the gravitomagnetic field then the tidal force — E"^X^ measured by the observers 
(i.e., the relative acceleration of two neighboring test particles of 4- velocity connected by X°) 
vanishes. If X" is orthogonal to the gravitomagnetic field, H"Xa = 0, then it is an eigenvector of 
Eq,^, of eigenvalue H^. The eingenvectors X" thus span a two dimensional subspace on the rest 
space ti* = 0, meaning that in these directions the tidal force is proportional to the separation 
vector X". Next we will physically interpret this for the special case of the Godel universe. 

4.1 The Godel Universe 

The Godel universe is a solution corresponding to an homogeneous rotating dust with negative 
cosmological constant. The homogeneity requires that the dust rotates around every point. The 
line element is given by: 

A^dx' = e^-^^dy , ^ijdx'dx^ = dx^ + ie^^'^^ciy^ + dz^ , (105) 

where w is a constant. It is straightforward to show that the gravitomagnetic field H is uniform, and 
equal to H = 2ujez; hence, by virtue of (jl03p . the gravitomagnetic tidal tensor vanishes: Ma/3 = 0. 
For this reason, this universe has been interpreted in [H |2] as being analogous to an uniform 
magnetic field in the curved 3-manifold with metric 7, and the homogeneous rotation physically 
interpreted in analogy with a gas of charged particles subject to an uniform magnetic field — in 
that case one equally has Larmor orbits around any point. 

Now we will interpret its gravitoelectric tidal tensor. In the coordinate system (I105|) it reads, 
for the = observers: 

It vanishes along z, and is isotropic in the spatial directions x, y orthogonal to H. It it similar 
to the Newtonian tidal tensor ViVjV of a potential: V = uJ^ix^ + y'^)/2, corresponding to a 2-D 
harmonic oscillator, which is the potential of the Newtonian analogue of the Godel Universe |91j : 
an infinite cylinder of dust rotating rigidly with angular velocity w. The potential V is such that 
the centrifugal force on each fluid element of the rigidly rotating cylinder exactly balances the 
gravitational attraction. This causes a curious effect in the Newtonian system: the fluid is seen 
to be rotating about any point at rest in the frame comoving with the "original" cylinder; indeed. 
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through an arbitrary point passes an axis of rotation relative to which the system is indistinguishable 
from the "original one". 

Therefore, whilst the gravitomagnetic field and tidal tensor, as well as the mapping via Klein- 
Gordon equation below, link to the magnetic analogue of the Godel universe, the gravitoelectric tidal 
tensor links to the Newtonian analogue, both yielding consistent models to picture the homogeneous 
rotation of this universe. 

5 Linear gravitoelectromagnetism 

The oldest and best known gravito-electromagnetic analogies are the ones based on linearized 
gravity, which have been worked out by many authors throughout the years, see e.g. [231 ESI [Ml 
[971 [Ml [211 [Zni 11061 l35l [7]. In the way it is more usually presented, one considers a metric given 
by small perturbations \hai3\ ^ 1 around Minkowski spacetime: gai3 = r]ai3 + ha^, and from the 
components /iq,/3 one defines the 3- vectors G and H, in terms of which the equations dubbed the 
metric is written in. The metric 

ds^ = -{1 + 2^>) dt"^ + 2Ajdtdx^ + [5ij + 2S,ij] dx'dx^ , (106) 

If ones considers stationary perturbations, as is more usual (e.g. [23l 1701 [T06l [351I7]). the GEM 
fields are (up to numerical factors in the different definitions): G = — V$, H = 'V x A, where 
Vj = d/dxi denotes herein (and only herein!) the spatial derivative operator associated to the 
background euclidean metric 5ij. These fields are straightforwardly related to the ones in Sec. [3l 
they are just, to linear order, minus the acceleration and twice the vorticity of the zero 3- velocity 
observers (u* = 0) with respect to the coordinate system associated to the metric (I106P (they can 
be called "static observers"). Thus they are simply a linear approximation to the Quasi-Maxwell 
fields in Eqs. ([50]). Note that, to hnear order, V$ V$, V x V x A 

If the fields depend on time, different definitions of the fields exist in the literature, as a complete, 
one to one GEM analogy based on inertial fields, holding simultaneously for geodesies and for the 
field equations, is not possible, as we shall see below (cf. also [981 [971 [Ml [21 [25]) • So if one chooses 
to write one of them in an electromagnetic like form, the other will contain extra terms. We stick 
to defining G and H by minus the acceleration and twice the vorticity of the congruence (i.e. the 
same definitions given in Sec. 13.21 for congruence adapted frames, only this time linearized), which 
seems to make more sense physically, as with this definition the fields appear in the equation of 
geodesies playing formally analogous roles to the electric and magnetic fields in the Lorentz force. 
That amounts to define: 

dA 

G = -V$-— ; H = VxA 
ot 

The space part of the linearized equation for the geodesies, in the coordinate basis ej, is obtained 
from the corresponding exact equation in the tetrad field ([45l) of Sec. [3] as follows. Let e"^^ denote 
the transformation matrix between the coordinate basis Bq, = d/da associated to the metric system 
in p06|) and the orthonormal tetrad adapted to the = observers; that is, = e^^Ba. To linear 
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order, e"^, and its inverse e"^ , are given by: 



- ^0 _ c i A 



eo = v^l + 2(/>eg; = ej + ^.^e-. + ^^e^ 

Thus, [/* = e*^ = f/* + iiUj\ putting also C/* = dx^ /dr, substituting in (j45]) . linearizing in the 
perturbations and in and noting that, to linear order, 

1 d^- dP 

= ^(^;,) = + « ^ = = ^ (107) 

the equation for the geodesies reads : 

-f2=G + ^;xF-2^i;^el (108) 

That is, the extra term compared to the Lorentz force in electromagnetism comes from the time 
derivative of the spatial metric (which is true also in the exact case, as we have seen in Sec. I3.2p . 
Noting that d'^x/dr'^ ~ d'^x/dt^ — vd^/dt, with v = dx/dt; we can also write this result as: 

l£ = o + 24.^5 + (100) 



The gravitational field equations in this regime are obtained by linearizing (I73p - (j78p and sub- 
stituting relations (|107p : 

V-G = -47r(2p + n)-^; (i) V x G = -— (ii) 

V-H = 0; (iii) V X i? = - 167rJ + 4^^.' ^'''^Sk (iv) (110) 

ot ■' 

Eqs. (IllOL ). (jllOL v). and (jllOl v). are, respectively, the time-time, time-space, and space-space com- 
ponents of Einstein's equations with sources (114k ): Eqs. (IllOt ii) and (IllOp are, respectively the 
time-time and space-time components of the identities (jl4b ). To obtain (jllOl v) from the exact 
Eq. (f72]) . we note that Rij reads, to linear order 

As for the time-space component of the identity (jl4b ). i.e., Eq. ()78p . it yields, to linear order, the 
trivial equation ^-R'^jj = 0. 

Eqs. (jllOp encompass two particularly important regimes: the "GEM limit", and gravitational 
radiation. Starting by the latter, in a source free region (T"^ = 0) one can, as is well known, 
through gauge transformations (employing the harmonic gauge condition, and further specializing 
to the transverse traceless, or radiation, gauge, see e.g. [TO]) make A = ^ = S,^^= C'' j = 0; with this 
choice, the only non trivial equation left is (jllO^ ^). yielding the 3-D wave equation d'^^ij/dt^ = V^^jj. 
The GEM regime is obtained making = (which effectively neglects radiation); in this case 
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the traceless shear of the congruence of zero 3-velocity observers = (in the coordmates system 
of (|106p ) vanishes: iTq^ = 0, and we have = 06ij/3 = —6ijd^/dt. This is the case also for the 
Post-Newtonian regime (e.g. \19 \ 147 1 I^O l 11161 169| ) ■ The somce is also assumed to be non-relativistic, 
so that the contribution of the pressure and stresses in Eq. (jllOp is negligible: 2p + ~ p. The 
equation for the geodesies (jlOSp then reads 



§ = G + v.H + 2'^v (111) 
at at 



and Eqs. (jllOp above become 



5^$ - dH 

V.G = -4^p + 3^; (i) VxG = -— (ii) 

BG r)^ A 

V-H = 0; (iii) V X # = - 167r/+ 4^1^ - 4 V? (iv) (112) 

ot ot^ 

a2f 



In some works, e.g. [24] . the gravitoelectric field is given a different definition: G' = -V^-\dA/dt; 
with this definition, and choosing the harmonic gauge condition, which implies V • ^ = —Ad^/dt, 
the non-Maxwellian term in Eq. (jll2t ) disappears; but, on the other hand, a "non-Lorentzian" 
term appears in the equations for the geodesies, where in the place of G in Eqs. (jl08p - (jl09p . we 
would have instead G' — ^dA/dt. The non-Maxwellian term in Eq. (|112[ v) is neglected in the Post 
Newtonian regime [471 [T9] . 

The presence of the terms dH/dt, "inducing" curls in G and H, respectively, analogous to 
the induction terms of electromagnetism, lead one to wonder if one can talk about gravitational 
induction effects in analogy with electrodynamics, and there is a debate concerning the applicability 
and physical content of this analogy for time-dependent fields, see e.g. [1] and references therein. 
Although a discussion of the approaches to this issue in the literature is not the scope of this work, 
still there are some points one can make based on the material herein. If one considers a time 
dependent gravitational field, such as the one generated by a moving point mass, e.g. Eq. (2.10) 
of |25j . indeed the gravitoelectric field G corresponding to this situation is different from the one 
for a point mass at rest, and has curl. That is, the acceleration —G of the congruence of observers 
at rest with respect to the background inertial frame (the "Post-Newtonian grid", e.g. [40]), gains 
a curl when the source moves with respect to that frame. From Eq. ()112t i). one can think about 
this curl as induced by the time- varying gravitomagnetic field H, see e.g. |116j . These fields 
are well suited to describe the apparent Newtonian and Coriolis-like accelerations of particles in. 
geodesic motion, as shown by Eq. (|llip — above (one must only bear in mind that in the case of 
time-dependent fields, the motion is not determined solely by G and H; there is an additional term 
with no analogue in the Lorentz force law, that leads to important significant differences). However 
the latter are mere artifacts of the reference frame; the physical (i.e., tidal) forces are a different 
story, one does not obtain the correct tidal forces by differentiation of these fields (as is the case 
with electrodynamics). Namely the curls of the GEM fields do not translate into these forces. The 
linearized gravito-electric tidal tensor, Eq. (195b ). reads in the GEM regime (-f^(jj) = —Sijd^/dt), 



1 dH'' 9$ 
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where we see that the curl (jll2t i) is subtracted from the derivative of G. That is, only the sym- 
metrized derivative G(j j) describes physical, covariant forces. This is manifest in that the curl of 
G does not induce a rotation on a set of neighboring particles (the gravitational field only shears 
them, see Sec. 12.21 and Eq. (f26|) therein), nor does it torque a rigid test body, see d]. Note that in 
electromagnetism, this rotation and torque are tidal manifestations of Faraday's Law of induction. 
Likewise, the curl of H is not manifest in the gravitomagnetic tidal effects (e.g., the force on a 
gyroscope); the linearized gravitomagnetic tidal tensor (195b ) reads, in this regime: 



1 

hi ~ — 



(114) 



where again we can see that the induction contribution 4dG/dt (and also the one of the term 
d'^A/dt^) to the curl of H is subtracted from the derivative of H. The physical consequences 
are explored in [3]: in electromagnetism, due to vacuum equation V x B = —dE/dt, there is a 
non-vanishing force on a magnetic dipole, F^j^^j = B^^fXj = V^{fl ■ B), whenever it moves in a non- 
homogeneous field (since the electric field it measures is time- varying, thus V x B ^ ^ Fem 7^ 0) . 
This expression for F^j^^ holds in the rest frame of the particle, for arbitrary fields. That is not 
necessarily the case in gravity. In vacuum, from Eqs. (I112t v) and (jll4p . we have Mij = — j)/2, 
and the gravitational force on a gyroscope, cf. Eq. ([ll2b) of Table [H is Fq = ^H^'^'^^Sj. Thus no 
analogous induction effect is manifest in the force, and in fact spinning particles in non-homogeneous 
gravitational fields can move along geodesies, which is exemplified in [1]. 

In the case that the field is stationary, we have a one to one correspondence with electromag- 
netism in inertial frames. Eq. (v) above becomes identical to (i) , and then we are left with a set 
of four equations — Eqs. (jllOn - (|110t v) with the time dependent terms dropped — similar, up to 
some factors, to the time-independent Maxwell equations in an inertial frame. These equations can 
also be obtained by linearization of Eqs. ([3}lb)-([3]4b) of Table [3l The space part of the equation 
of the geodesies: dU /dt = G + v x H , ci. Eq. (jlOSp above, is also similar to the Lorentz force in 
a Lorentz frame. The equation for the evolution of the spin vector of a gyroscope becomes simply 
dS/dr = S X H /2, which gives the precession relative to the background Minkowski frame, and is 
similar to the precession of a magnetic dipole in a magnetic field. The force on a gyroscope whose 
center of mass it at rest is Fq = 'V{S-H)/2, similar to the force Fem = V(/I--B) on a magnetic dipole 
at rest in a Lorentz frame; the same for the differential precession of gyroscopes/dipoles at rest, they 
read respectively, for a spatial separation vector Sx", 6Qg = — V(5x • H)/2, S^Iem = —V{6x- B). 



6 The formal analogy between gravitational tidal tensors and 
electromagnetic fields 

There is a set of analogies, based on exact expressions, relating the Maxwell tensor and the 
Weyl tensor Cafs^s, namely: 1) they both irreducibly decompose in two electric and magnetic type 
spatial tensors; these tensors obey differential equations — Maxwell's equations and the so called 
"higher order" gravitational field equations — which are formally analogous to a certain extent 
[30\ \3T\ [85| [32I [51], and they form, moreover, invariants in a similar fashion to the relativistic 
invariants formed by the electric and magnetic fields [30^ [31} \52\ [53] . In this section we will briefly 
review this analogy and clarify its physical content in the light of the previous approaches. 
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The Maxwell tensor splits, with respect to a unit time- like vector in its electric and magnetic 
parts: 

= F"^f//5 , = i.F'^pU^ , (115) 

i.e., the electric and magnetic fields as measured by the observers of 4- velocity U". These are spa- 
tial vectors: E'^Ua = B°'Ua = 0, thus possessing 3-1-3 independent components, which completely 
encode the characterize the 6 independent components of i^^i/, as can be seen explicitly in decom- 
positions ([I]) . In spite of their dependence on , one can use Ea and Ba to define two tensorial 
quantities which are independent, namely 

E'^E^-B-B^ = --^^ , £;"i?. = - ; (116) 

these are the only algebraically independent invariants one can define from the Maxwell tensor. 
The Weyl tensor has a formally similar decomposition: with respect to a unit time-like vector 
it splits irreducibly in its electric £ai3 and magnetic Hap parts: 

These two spatial tensors, both of which are symmetric and traceless (hence have 5 independent 
components each), completely encode the 10 independent components of the Weyl tensor, as can 
be seen by writing [32] 

C^/-^ = 4 {2C/[,[/[^ + g^;^] + 2 {e^p.^U^m'^'^U'' + e^^'^'^U^^n^^^U,} . (118) 

Again, in spite of their dependence on one can use £ai3 and Tia/B to define the two tensorial 
quantities which are f/" independent 

s^^s^p - n'^^nafs = ^^^^ — , £''^n^p = "^^'^^ ^ — . (ii9) 

8 16 
which are formally analogous to the electromagnetic scalar invariants ()116p . However, by contrast 
with the latter, these are not the only independent scalar invariants one can construct from Cap^v] 
there are also two cubic invariants, see [631 ED [Ml E5]. 

As stated above, these tensors obey also differential equations which have some formal similar- 
ities with Maxwell's; such equations, dubbed the "higher order field equations" are obtained from 
Bianchi identities RtjT[fj.u;a\ = 0- These, together with field equations (jl4b ). lead to: 

C^Va* = Svr (^T,[,.,] - i<7.[rr;,]) , (120) 

Expressing Caps^y in terms of £a(3, using (jllSp . and taking time and space projections of 
(jl20p using the projectors ([2|), we obtain, assuming a perfect fluid, the set of equations 

(121) 

Dp 
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where e^^p = ^puprU'^, V denotes the spatially projected covariant derivative, cf. Eq. ([56|) . curlAap = 
('^'^(^^Ap^y.p^, and the index notation (^z>) stands for the spatially projected, symmetric and trace 
free part of a rank two tensor (cf. definitions in |32j): 

^) (^iiu = Di^pUy), bj" = e°^p^U^-p/2, a" are, respectively, the expansion, shear, vorticity, and 
acceleration of the congruence of observers with 4- velocity C/". 

The analogous electromagnetic equations are the ones in Sec. 13.4. l) 

V^E^" = Anpc + 2ujpB'' (123) 
£-7/3^^.^ = _ + IffE'' - e^co^E^ + e%^B^a? + 47rj<"> (124) 

V^S^^ = -2ujpEi' (125) 

e'^^^Ep,^ = -^^^ + <y%BP - loB'' + e^t^^i?^ + e^^'^i^^a. (126) 

Eqs. p23p and p25p follow from Eqs. (I54p and ()58p . respectively, by noting that, for an arbitrary 
vector A"', 

A, = (t^ + hi) (t^ + hi) A\ = {Tl + hi) A]^ = A^a, + V^A^. 
Eqs. (fT2i]) and (fT26P follow from Eqs. ^ and ([Ml) by decomposing = 0-^/3 + Ohap/^. 
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6.1 Matte's equations vs Maxwell equations. Tidal tensor interpretation of 
gravitational radiation. 



Table 4: Formal analogy between Maxwell equations (differential equations for electromagnetic 
fields) and Matte's equations (differential equations for gravitational tidal tensors) 

Electromagnetism Linearized Gravity 



Maxwell's Equations 



Matte's Equations 








(Ula) 







(Hlb) 







(Il2a) 







(|l2b) 






gl3a) 




dm 


(|l3b) 


dt 


= - 


dt 




dt 


(Il4a) 




dt 


(jl4b) 



Wave equations 
52 



d^dk ]E' = dlSa) 



dt^ 



^ - d^du ) = (Il6a) 



Wave equations 

d'^dk ) E,, = 



(jl5b) 
gl6b) 



In vacuum, the Bianchi identities become: 

= 0; (a) «%-,;^ = (6) (127) 

(the second equation following from the first and from vacuum equation R^^^ = 0). The formal 
analogy with Eqs. ([6]), for = 0, is now more clear [31] • In a nearly Lorentz frame, and to 
linear order in the metric potentials, Eqs. (|12ip - (jl22p . for vacuum, become Eqs. ([311b)-(Ill4b) of 
Table m which are formally similar to Maxwell's equations in a Lorentz frame (Illla)-(jll4a). The 
analogy in Eqs. ([Ill)-(I114) was first found by Matte |30j . and further studied by some other 
authorfl [SH ED [62]. Taking curls of Eqs. ([I3a)-([l4a) we obtain the wave equations for the 
electromagnetic fields; and taking curls of ((H3a)-([3l4a), we obtain gravitational waves, as wave 
equations for gravitational tidal tensors. 

Hence, to this degree of accuracy, vacuum gravitational waves can be cast as a pair of oscillatory 
tidal tensors Eq,^, Hq,^ propagating in space by mutually inducing each other, just like the pair 
of fields E"-, in the case of the electromagnetic waves. Also, just like E°' and B" are equal 
in magnitude and mutually orthogonal for a purely radiative field, the same applies to the waves 
in ((4l5b)-((4j6b) of Table [H In the electromagnetic case this implies that the two invariants (1116P 
vanish; likewise, the gravitational invariants (jllOp also vanish for a solution corresponding to pure 

*The exact equations (|12ip - (|122p . in vacuum, can take the form (|3]lb)-(|3]4b) of Tabled in a local Lorentz frame, 
where a" = uj°' — Oap = = 0, as done in e.g. [102) . That however holds only in a small neighborhood around a 
particular observer worldline, and is therefore not suited to study the propagation of radiation. 
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gravitational radiation according to Bel's second criterion (cf. e.g. [60] p. 53) — a definition based 
on "super-energy", discussed below. 

An interesting aspect of this formulation of gravitational radiation, contrasting with the more 
usual approaches in the literature, e.g. [TOl [28l [231 [3] — which consist of equations for the prop- 
agation of gauge fields (the components of the metric tensor) having no local physical significance 
(only their second derivatives may be related to physically measurable quantities, see in this re- 
spect [75]) — is that Eqs. (|4l5b)-([4l6b) are equations for the propagation of tensors of physical 
forces, with direct translation in physical effects - the relative acceleration of two neighboring 
test particles via geodesic deviation equation ([Hlb), and the force on a spinning test particle, via 
Mathisson-Papapetrou-Pirani Eq. ([Tl2b), or the relative precession of two nearby gyroscopes, via 
Eq. (j27p . The latter effects come from the "magnetic part" (i.e., involving Hq,^) of the radiation, 
little-studied in comparison to its electric counterpart; its role, in equal footing with the electric 
part, is made explicit in this formulation — they are both essential for the existence of gravitational 
waves, as pair of mutually inducing fields is needed. It yields a simple and intuitive description 
of the interaction of a gravitational wave with a pole-dipole spinning particle, as a coupling of the 
wave equation ([4l6b) to the spin- vector of the particle, via Eq. ([Tl2b). In other words: putting 
together Matte's equations with the physical interpretation of the electric and magnetic parts of 
the Riemann tensor given in Sec. [21 we obtain this suggestive interpretation of gravitational ra- 
diation, which has not (to our knowledge) been put forth in the existing literature. That might 
be down to two reasons: on the one hand the works treating the higher order/Matte equations 
[30l[3ll[Z5l[85l[6ll[62l[5ll[87],as well the tensors {£ai3,'Hai3} m other contexts e.g. [8l[8Zl[56], lack 
a physical interpretation of the tensors Tiai3 /^ai3 (the tensors £ai3 /^ai3 in their turn are reasonably 
well understood due to their role in the geodesic deviation equation), which are portrayed as not 
well understood [84[ [871 [56] , or given inconsistent interpretations that immediately lead to contra- 
dictions [82 [ I83 [ [32] (see [2] for more details). On the other hand, in the literature concerning the 
interaction of gravitational waves with spinning particles, e.g. [72 [l71l[73l l74j. neither gravitational 
waves are written in the form ([4l5b)-([4l6b), nor, as a matter of fact, is the force on a spinning 
particle explicitly related to the magnetic part of the Riemann tensor Mai3 (or to Tiajs), which is 
likely down to the fact that in these treatments the Tulczyjew-Dixon spin condition 3°^^ Pj^ = is 
employed in the Mathisson-Papapetrou equation (first term of Eq. ()130p ). instead of the Mathisson- 
Pirani condition 3°^^ Up = 0. The two conditions are equally valid (see in this respect [5]); however 
it is only when one uses the latter, that one obtains the expression ([Tl2b) for the force (see [1] for 
more details). 

It is also interesting to note that in the traditional treatments the wave equations are obtained 
from a linearized form of Einstein equations ()14h ): whereas equations ((3l5b)-([3l6b) come from a 
linearization of the higher order field equations (jl27b ). (Even though the former still play a role, 
as in order to obtain Eqs. ([4l5b)-([4l6b) from the differential Bianchi identity (jl27h ). we have 
substitute Rai3 by the source terms using Einstein's equations ([T4k )). 

It is important to realize that whereas in the electromagnetic radiation are the vector fields that 
propagate, gravitational radiation is a purely tidal effect, i.e., traveling tidal tensors not subsidiary 
to any associated (electromagnetic-like, or Newtonian-like) vector "field"; it is well know that there 
are no vector waves in gravity (see e.g. [651 [28] I64| : such waves would carry negative energy if 
they were to exist, cf. [28] p. 179). We have seen in Sec. 13.51 except for the very special cases of 
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the linear regime in weak, stationary fields (and static observers therein), the gravitational tidal 
tensors cannot be cast as derivatives of some vector field. In the electromagnetic case of course there 
are also tidal effects associated to the wave; but their dynamics follows triviall}{§ from Eqs.(jll3a)- 
(jll4a); to this accuracy, the tidal tensors as measured by the background static observers are just 
Eij = Eij, Bij = Bij; hence the equations of their evolution (i.e., the "electromagnetic higher 
order equations") are: 

= e>'S,7,fc = (128) 

^ij,k — ^ik,j — — ^ij,k — -D/fcj — —g^ U^yj 

These four equations are the physical analogues of the pair of gravitational Eqs. ([Il3b)-(4b); we 
have two more equations in electromagnetism since Eij and By are not symmetric. Eqs. (jl28p . 
and the first equality in Eqs. (jl29p . come from the fact that derivatives in flat spacetime commute; 
therefore e^^^Eji^k = (-i^^Ej^yij^-^ = and Eij^k = -E'zjfe = ^ikj- Thus, Eqs. (jl29p . which are the only 
ones that contain dynamics, are obtained by simply differentiating d/dj Eqs. (j3l3a)-(j3J4a). And 
the wave equations for the electromagnetic tidal tensors follow likewise from differentiating d/dj 
Eqs. (j45a)-([4t)a). Note also that the fact that in gravity, IHIj[«,fc] / is again related to the fact that, 
even in the linear regime, the gravitational tidal tensors are not derivatives of some vector fields. 

The tidal tensor interpretation of gravitational waves gives insight into some of their fundamen- 
tal aspects; we will mention three. Firstly, being gravitational waves (in vacuum) traveling tidal 
tensors, then according to discussion in [21 they couple to dipole particles (i.e., spinning particles, or 
"gravitomagnetic dipoles") causing a force, and they can also cause a relative acceleration between 
two neighboring monopole particles; but they cannot exert any force on a monopole test particle 
(by contrast with their electromagnetic counterparts); that is the reason why a gravitational wave 
distorts a (approximately monopole) test body, but does not displace its mass center, as is well 
known, e.g. |23j . Let us state this important point in terms of rigorous equations. As explained in 
Sec. O an extended test body may be represented by its electromagnetic and "gravitational" mul- 
tipole moments (i.e., the moments of and Tp^ , respectively). The force exerted on a charged 
body in a electromagnetic field is, up to quadrupole order [9l [35| H]: 

^^^^ = qF-^U^ + \f^^-^^Q,, + \Qp,sF^''^- + ... 

where is the particle's 4-velocity, q is the charge, Qap is the dipole moment 2-form, which 
we may write as Qq,^ = Id^^Uj^T^ + ea/B-ySfJ-'^U^ , and Qa/s^y is a quadrupole moment of the charge 
4-current density Using decomposition ([7|), we can write the second term in terms of tidal 
tensors {F^^'^''°^Q^^/2 = El^'^dp + B^'^^p), the third term in terms of derivatives of the tidal tensors, 
and so on. Thus we see that in the electromagnetic case, there are force terms coming from the 
electromagnetic monopole, dipole, quadrupole moments, etc, corresponding to the coupling to, 
respectively, the field, the tidal field, derivatives of tidal field, etc. 

The exact equations of motion for a test body in a gravitational field are [8l[35], up to quadrupole 
order: 

= --R-^^^S'^'^U^ - -J^.s.R^^'''-''' + ... (130) 



We thank J. Penedones for discussions on this point. 
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where Sap = ^a/B^ysS^U^ is the spin 2-form, which is the dipole moment of the mass current 4- vector 
Jp = —Tp^Up, see e.g. [311], and Jp^Sa is a quadrupole moment of the energy- momentum tensor 
of the test particle, Tp^ . The first term yields of course — HI'^"S'/3, in agreement with Eq. ([Tj2b), 
corresponding to the coupling of the gravitomagnetic tidal tensor with the spin vector S"" (the 
gravitomagnetic dipole moment), the second term the coupling of the derivatives of the tidal field 
to the quadrupole moment, and so on. Hence we see that, by contrast with electromagnetism, in the 
gravitational case the interaction starts only at dipole order; as the lowest order (in differentiation) 
physical fields present in gravity are the tidal tensors (the real, physical content of the gravitational 
waves), which do not couple to monopoles. 

Secondly, it sheds also some light on the issues of gravitational induction. As seen in Sees. [21 
13.51 and [5l electromagnetic-like induction effects are absent in the gravitational physical forces. 
However, instead of some vector fields, in gravity one can talk of a different type of induction 
phenomena, for tidal tensors: a time- varying magnetic tidal tensor induces an electric tidal tensor, 
and vice- versa, as implied by Eqs. ([Il3b)-([3l4b). These tidal tensors propagate by mutually inducing 
each other, giving rise to the gravitational radiation, just like the laws of electromagnetic induction 
lead to electromagnetic radiation. 

Finally, one can also relate the fact that in the gravitational waves what propagates are the 
tidal tensors, without any associated electromagnetic like vector field, with the fact that, to the 
gravitational waves (and the gravitational field itself), one cannot associate a local energy or mo- 
mentum density, in the traditional sense of being manifest in the energy momentum tensor T^P . 
This becomes more clear contrasting with the electromagnetic case; the electromagnetic field gives 
a contribution to the total energy momentum tensor T"^ equal to: 



t; 



«/3 



EM 



1 r 
8^ 



(131) 



leading to an energy density pem and spatial momentum density (Poynting vector), with 

respect to an observer of 4- velocity u": 

PEM = T'^iruaup = ^[{E^r{E-)^ + {B-rma\ (132) 
A = -T^^l^up = ^e%,^u%E-r{B^r (133) 

or, in vector notation, SnpEM = E{u)'^ + B^u)"^; pem = E{u) x B{u)/4tt. Now note that the 
significance of these expressions as energy and momentum densities can be traced back |76[ [26] , 
at the most fundamental level, to the Lorentz force (together with Maxwell equations), the work 
done by it, and its spatial momentum transfer. It is thus clear that it cannot have a direct physical 
gravitational analogue, as there is no physical, covariant, gravitational counterpart to the Lorentz 
force and the vector fields E, B; the inertial "force" MSh . and its Lorentz- like form for stationary 
spacetimes ([48]) (both being but the geodesic equation in a different language, no real force being 
involved), as well as the fields G and H, are mere artifacts of the reference frame, which can be 
made to vanish by switching to a locally inertial frame, as explained in Sec. [3l It is the same 
for the "energy" and "momentum densities" arising from the Landau-Lifshitz pseudo-tensor t^j,, 
e.g. [77\ 128]: the "momentum density " arising from t^^, to Post-Newtonian order, may actually be 
written as [35]: ~ {-G x H + 3G6)/47r 
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6.1.1 Super-energy 

On the other hand there is a quantity built on tidal tensors (not on vector fields), having thereby 
a local physical existence, that seems to fit well in the tidal tensor interpretation of gravitational 
radiation — the so-called Bel-Robinson super-energy tensor [31], which reads in vacuum: 



2 L 



(134) 



(the more general expression in the presence of sources is given in e.g. [78]; a general superenergy 
tensor can also be defined, e.g. [321 [78] from the Weyl tensor; it is obtained by simply replacing 
RaP'yS — ^ Rai3"/s)- The formal analogy to the electromagnetic energy- momentum tensor, Eq. (I13ip 
above, is clear. The super-energy density and flux are defined as: 

W ^ T'^P^'u^upu^us = ^ [(E«)-^(E"),;3 + (IH")"^(IHI")„;3] (135) 
p(a) _ -T^o.)p^5^^^^^^ ^ e''^^^u''{WY^{W)\ (136) 

where the gravitational tidal tensors play a role formally analogous to the electromagnetic fields 
in the quantities pem and Eqs. (|132p - (|133|) above. Note that expressions (|135p - (|136|) are 

exact. W is positive definite, and is zero only if the curvature vanishes [311120]; it thus qualitatively 
agrees with one would expect from field energy associated with the physical gravitational field (i.e., 
the curvature tensor). According to Bel |3H I60|. gravitational radiation is characterized by a flux 
of super-energy "P" (which, at least in the linear regime, is clear from the equations ([lllb)-([3l6b) 
above), parallel to the spatial direction of propagation of the wave. And that, at a given point, the 
non-vanishing oiV^"^ for every observer tt" (i.e., an intrinsic super-flux), is sufficient to ensure that 
gravitational radiation is present. Several criteri^^ have been proposed (see jllOj for a review) to 
characterize radiative states, based on the electromagnetic analogy. It is indeed tempting to think 
of super-energy and super-momentum (or super-flux) as the form in which gravitational waves carry 
what then, in the interaction with matter — e.g., through the test body multipole moments, as 
described by Eq. ()130p — manifests itself as ordinary energy and momentum. This is an hypothesis 
to be studied in detail elsewhere; herein we would just like to point out that this is line with the 
point of view in [75] (where it is also shown that, via W , gravitational radiation has an active 
attractive effect), and that the fact that super-energy is not a conserved quantity in general in the 
presence of matter [79[ I78j. seems also to strenghten the hypothesis. There are however conceptual 
difficulties [78] [HTl I109j in the physical interpretation of these quantities, such as the "strange" 
dimensions of W , which are that (taking L = M) can be interpreted an energy density per 
unit area [Ml [TH [79] , or of an energy density times a frequency squared [75] , or an energy density 
squared [78]. The very question if the superenergy has any physical reality is an open question, as is 
still the old problem of the definition of the energy of the gravitational field (which we will not try 
to further elaborate herein). In the spirit of the arguments given above — that the GEM fields are 
frame artifacts, the non-existence of local quantities physically analogous to the electromagnetic 
vector fields (consequence of the equivalence principle), and that gravitational waves are tidal 
tensors propagating by mutually inducing each other (without associated vector fields), which do 
not interact with the multipole structure of test bodies in the same way electromagnetic waves do 



^°We thank A. Garcia-Parrado for discussions on this issue. 
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(in particular, do not couple to monopole particles) — it seems plausible that the dimensions of 
the fundamental "entity" carried by gravitational waves are not the same of their electromagnetic 
counterparts, hence not the ones of an energy. This is consistent with the point of view in |49 1ll09j . 

6.2 The relationship with the other GEM analogies 

The analogy drawn in this section is between the electromagnetic fields and the electric and magnetic 
parts of the Weyl tensor: {E"', S"} -f-)- {£^'^ , H^"^. It is clear, from the discussion of the physical 
meaning of {E^'', H'"^} in Sec. [21 and from the dynamical gravitational counterparts of {E°^, S"} 
discussed in Sec. [3l that this analogy is a purely formal one. It draws a parallelism between 
electromagnetic fields, whose dynamical gravitational analogues are the GEM inertial fields {G, H} 
of Sec. O with gravitational tidal fields, which are the physical analogues not of {E°', B"'}, but 
of the electromagnetic tidal tensors {-Ea^, -Ba^} (which, for an observer of fixed 4- velocity U", 
are covariant derivatives of the E" and it measures). This sheds light on some conceptual 
difficulties in the literature regarding the physical content of the analogy and in particular the 
physical interpretation of the tensor H^'^ , see [2j for details. And is also of crucial importance for 
the correct understanding of physical meaning of the curvature invariants, and their implications 
on the motion of test particles, which is subject of detailed study in [63]. 

7 When can gravity be similar to electromagnetism? 

A crucial point to realize is that the two exact physical gravito-electromagnetic analogies — the 
tidal tensor analogy of Sec. [2| and the inertial GEM fields analogy of Sec. [2] — do not rely on a 
close physical similarity between the interactions; as the gravitational objects {G, E^^, Hq^}, 
despite playing analogous dynamical roles to the objects {E, B, E^p, -Bq,^} in electromagnetism, are 
themselves in general very different from the latter. Even for seemingly analogous setups (e.g. the 
EM field of spinning charge, and the gravitational field of a spinning mass). In this sense these 
analogies have a different status compared to the popular GEM analogy based on linearized theory, 
which, in order to hold, indeed requires a close similarity between the two interactions, to which 
the former two are not bound. 

What the tidal tensor formalism of Sec. [21 together with the inertial fields formalism of Sec. [3l 
provide, is a "set of tools" to determine under which precise conditions a similarity between the 
gravitational and electromagnetic interactions may be expected. 

The key differences between electromagnetic and gravitational tidal tensors are: a) they do 
not exhibit, generically, the same symmetries; b) gravitational tidal tensors are spatial whereas 
electromagnetic ones are not; c) electromagnetic tidal tensors are linear, whereas the gravitational 
ones are not. 

The electromagnetic tidal tensors, for a given observer, only have the same symmetries and 
time-projections as the gravitational ones when the Maxwell tensor is covariantly constant along 
the observer's worldline; that is implied by Eqs. ((Il4) and ([Il6). That restricts the eligible setups 
to intrinsically stationary (i.e., whose time-dependence, if it exists, can be gauged away by a 
change of frame) fields, and to a special class of observers therein; for electromagnetic fields in 
flat spacetime, those observers must be static in the inertial frame where the flelds are explicitly 
time-independent. This is an important point that is worth discussing with some detail. Consider 
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the two basic analogous fields, Coulomb field of a point charge, and Schwarzschild gravitational 
field. Consider also in the latter observers O in circular motion: 4-velocity C/" = {U^ ,0,0,11'^), 
angular velocity Q = U'^/U'^; the worldlines of these observers are tangent to killing vectors: 
[/" II C^QajB = 0; hence one uses to say (e.g. [281 EZ]) that they see a constant spacetime 
geometry; for this reason they are called "stationary observers". Now consider observers in circular 
motion around a Coulomb charge. Despite moving along worldlines tangent to vector which is 
a symmetry of the electromagnetic field: C/" || C^F^fi = 0, the observers C/" do not see a 
covariantly constant field: Fap-^V^ ^ (they see a constant magnetic field and an electric field 
E constant in magnitude but varying in direction), which by virtue of Eqs. (HJGa), ([Tj4a), implies 
that the electromagnetic tidal tensors have an antisymmetric part (in particular the spatial part 
^[ij] 7^ 0)> ^'^d thus means that they cannot be similar to their gravitational counterparts. This 
is a natural consequence of Maxwell's equations and can be easily understood as follows. The 
magnetic tidal tensor measured by O is a covariant derivative of the magnetic field as measured in 
the inertial frame momentarily comoving with O: Bap = -kFary-f^U"^ = Ba-i3\u=const. = {BMCRF)a;i3- 
Now, -Bjjj] 7^ means that Bmcrf has a curl; which is natural since in the MCRF the electric field 
is time-dependent, which, by virtue of Maxwell equation V x B = dE/dt = ^E x Q, (holding in the 
MCRF, and for which ([Tj6a) is a covariant form) induces a curl in B. 

But even if one considers static observers in stationary fields, so that the GR and EM tidal 
tensors have the same symmetries, still it does not mean a close similarity between the interactions. 
The electromagnetic tidal tensors are essentially derivatives of the electromagnetic fields, being 
themselves differential fields linear in the electromagnetic 4-potential A" = ((/>, A), whereas the 
gravitational ones are non-linear in the GEM fields, as shown by Eqs. (j87p - (j88p . being the grav- 
itomagnetic field H itself non-linear in the metric potentials A. This means that one expects 
a similarity between tidal tensors in the two limiting cases — linearized theory, and the ultrasta- 
tionary spacetimes considered in Sec. llOOl where $ = G = 0, and, therefore, cf. Eqs. ([88j) and ([50|) . 
the exact gravitomagnetic tidal tensor is linear: H--- = -D-(V x A)i/2. We have seen in Sec. llOOl 
that there is indeed an exact mapping (via Klein-Gordon equation) between the dynamics in these 
spacetimes and an electromagnetic setup. 

To what concerns the concrete effects, the precise conditions (namely regarding the time de- 
pendence of the fields) for occurrence of a gravito-electromagnetic similarity are specific to the 
type of effect (see also [25]). For the covariant effects (implying physical gravitational forces, i.e., 
tidal forces) such as the force on a spinning particle gyroscope or the worldline deviation of two 
neighboring particles, are the tidal tensors as measured by the test particles of 4-velocity that 
determine the effects, cf. Eqs. ([Ill)-([Il2); which means that it is along the particle's worldline that 
the time independence is required. That basically implies that the similarity only occurs in the 
instant when the particles are at rest in stationary fields, so it does not hold in a truly dynamical 
situation. For the case of the correspondence between the Lorentz force Eq. (j46p . and the geodesic 
equation formulated as an inertial force, which is a coordinate effect, from Eq. ()45p we see that 
the requirement is that the frame is rigid, i.e. ctq,^ = ^ = 0; as explained in Sec. 13.21 that amounts 
to say that the spatial part of the metric (in the coordinates associated to such frame) must be 
time-independent. This can also be stated in the following manner, generalizing to the exact case 
the conclusion obtained in [25] (in the context of the Post-Newtonian approximation): in the case 
of the GEM analogy for the geodesic equation, the stationarity of the fields is required in the ob- 
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server's frame (not in test particle's frame! The test particles can move along arbitrary worldlines). 
As for the gyroscope "precession" (j52p and the correspondence with the precession of a magnetic 
dipole (j53p . there is no restriction on the time dependence of the fields. 

8 Conclusion 

In this work we collected and further developed different gravito-electromagnetic analogies existing 
in the literature, and clarified the connection between them. A detailed summary of the material 
in this paper is given in the introduction; we conclude by briefly summarizing the main outcome 
of each approach, and their applicability. The analogies split in two classes: physical and purely 
formal. In the second category falls the analogy between the electric and magnetic parts of the 
Weyl and Maxwell tensors, discussed in Sec[6l The physical analogies divide in two classes: exact 
analogies, and the best known Post-Newtonian and linearized theory approaches. Exact physical 
analogies are the analogy between the electromagnetic fields and the inertial fields of Sec. [3l and 
the tidal tensor analogy of Sec. [21 

These analogies are useful from a practical point of view, as they provide a familiar formalism, 
and insight from electromagnetic phenomena to describe otherwise more complicated gravitational 
problems. Indeed, there is a number of fundamental equations, summarized in Tabled which can 
be obtained from the electromagnetic counterparts by simple application of the analogy. But the 
existence of these analogies, especially the exact, physical ones, is also interesting from the the- 
oretical point of view, with intriguing similarities — both in the tidal tensor, and in the inertial 
field formalism, manifest in Tables [T] and O respectively — and enlightening differences unveiled 
(namely the ones manifest in the symmetries of the tidal tensors). The deep connection between 
gravitation and electrodynamics is still very much an open question; the similarities and fundamen- 
tal differences unveiled in the analogies and formalism herein may be a small step in the direction 
of that understanding. 

The tidal tensor formalism is primarily suited for a transparent comparison between the two 
interactions, since it is based on mathematical objects describing covariant physical forces common 
to both theories. The most natural application of the analogy is the dynamics of spinning multipole 
test particles, which is studied in the companion paper [3] to quadrupole order. As a formalism, 
consisting of mathematical objects with direct physical interpretation, encoding the gravitational 
physical, covariant forces, it can be useful in many applications, namely gravitational radiation (as 
discussed in Sec. l6.ip . and whenever one wishes to study the physical aspects of spacetime curvature 
(as a further example, we mention the physical significance of the curvature invariants, which have 
attracted attention in the context of experimental astrophysics, and are studied in this framework 
in [63]). 

The analogy based on inertial GEM fields from the 1+3 formalism, Sec. El is a very powerful 
formalism, with vast applications; especially in the case of stationary spacetimes, where for arbi- 
trarily strong fields the equation for geodesies is cast in a form similar to Lorentz force; many other 
effects related to frame dragging can be treated exactly with the GEM fields: gyroscope "precession" 
(relative to the "distant stars") d (H [201 iQl E] , the gravitomagnetic clock effect [PUT REF], the 
Sagnac effect |112j . the Faraday rotation [9l], the force on a gyroscope (Sec. 13.61 and [H]; note 
however that it is not as general as the the tidal tensor formulation of the same force); and other 
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applications, such as the matching of stationary sohitions [16], or describing the hidden momentum 
of spinning particles [1]. The general formulation of GEM fields in Sec. [3l applying to arbitrary 
fields and frames, extends the realm of applicability of this formalism. 



Table 5: What can be computed by direct application of the GEM analogies 



Result 


Approach 


• Geodesic deviation equation ([Hlb): 

"Til" fTT'lif TTT'I" /fn 1 \ 

-Replacmg {q,hai3\ — )■ {m, — Ito/?) m dUlaj. 

• Force on a gyroscope ([TJlb): 

-Replacing {/x",^^/?} ^ {S", -Mo,^} in ([Ula). 

• Differential Precession ([Tj8b): 
-Replacing {a,Ba/3} — )• {!,— Hq,^} in ([Ula). 

• Temporal part of Einstein's equations ([HSb) -([Tj6b): 
-Replacing {Ea,3,B^fs} ^ {E„/3,E[„^} in Eqs. 

and pc^2p + in 1^, j° 2J° in (fT2]). 


Tidal tensor analogy 

(Exact, general results) 


• Geodesic Equation (|18]) {stationary fields) 
-Replacing {E, B} — )• {G, H} in (06]), multiplying by 7. 

• Gyroscope "precession" Eq. ([52]) {arbitrary fields): 
-Replacing {p,B} {S,H/2} in ([53]). 

• Force on gyroscope Eq. (I98p {stationary fields, 
particle's worldline tangent to time-like Killing vector): 

-Replacing {fl,E,B} {S,G,H/2} in ([96]), factor 
of 2 in the last term. 


Inertial "GEM fields" analogy 

(Exact results, require special frames) 


• Higher order field equations ([4llb)-([4j4b): 
-Replacing {E,B} {Eij,]Hij} in Eqs. ([3la)-([l4a). 

• Equations of gravitational waves ([Il5b)-([3l6b): 
-Replacing {E,B} {Ey,EIjj} in Eqs. ([I5a)-([l6a). 


Weyl- Maxwell tensors analogy 

(Results for linearized theory) 



The well known analogies between electromagnetism and Post-Newtonian, and linearized grav- 
ity, follow as a limiting case of the exact approach in Sec. [S] In the case of the tidal effects, they 
can be seen also as limiting case of the tidal tensor analogy of Sec. [2] (in the sense that for weak, 
time-independent fields, the tidal tensors are derivatives of the GEM fields). Acknowledging this 
fact, and understanding the conditions under which linear GEM is obtained from the rigorous, 
exact approaches, is important for a correct interpretation of the physical meaning of the quantities 
involved, which is not clear in the usual derivations in the literature (this is especially the case for 
many works on linear GEM), and thus prone to misconceptions |2l [1]. On the other hand, linear 
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GEM is the most important in the context of experimental physics, as it pertains all gravitomag- 
netic effects detected to date [HI 11061 HH 11161 11151 1114j , and the ones we hope to detect in the 
near future [13] , 

As for the analogy between the electric and magnetic parts of the Weyl and Maxwell tensors, its 
most important application is gravitational radiation, where it provides equations for the propaga- 
tion of tensors of physical forces (not components of the metric tensor, as in the usual approaches, 
which are pure gauge fields) , with direct translation in physical effects via the tidal tensor formalism 
of Sec. [21 This analogy has also been used to address the fundamental questions of the content 
of gravitational waves, and the energy of the gravitational field. Namely, to define covariant, lo- 
cal quantities alternative to the gravitational energy and momentum given by the Landau-Lifshitz 
pseudo-tensor (which can only have a meaning in a global sense, and in asymptotically flat space- 
times): the super-energy and super-momentum encoded in the Bel tensor. The motivation for the 
definition of this tensor is the analogy with electromagnetism; and the existing criteria for radiative 
states |110| . states of intrinsic radiation [311 I109j or pure radiation ( |113| . see also [60] p. 53), 
are also solely driven by it. The analogy is also useful for the understanding of the quadratic 
invariants of the curvature tensor; indeed, it will be shown elsewhere [63] that using the two ap- 
proaches together — the formal analogy between the Weyl and Maxwell tensors to gain insight 
into the invariant structure, and the tidal tensor analogy as a physical guiding principle — one can 
explain, in the astrophysical applications of present experimental interest (as mentioned above), 
the significance of the curvature invariants and the implications on the motion of test particles. 
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